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1. Introduction 

The purpose of these lectures is to present a short review of some recent developments 
of the lattice approach to the problem of the quantization of gravity. This approach 
goes back to the seminal paper by Regge (Regge (1961)), but many new and sometimes 
surprising results have been obtained during the last ten years, thanks both to the 
progresses of computers and of numerical techniques, and to new theoretical ideas, 
developed in particular for two dimensional (2d) theories of gravity and for string 
theories. These developments came also from a convergence of different fields of 
theoretical physics: quantum gravity, string and superstring theories, but also lattice 
gauge theories and statistical mechanics, in particular random surface models. 

The construction of a theory of gravity which encompasses both General Relativity 
(GR), or some of its extensions, and the principles of Quantum Mechanics is probably 
one of the most difficult and exciting problems of theoretical physics. It is well-known 
that there are at least two basic difficulties. Firstly, the coupling constant for GR 
(Newton's constant) has the dimension of the inverse of a mass squared, and leads 
to a non-renormalizable quantum theory, with violations of unitarity at Planckian 
energies. Secondly, the Einstein equations (the evolution equations for the metric) 
are derived from an action principle, but the Einstein-Hilbert action is unbounded 
from below, and a naive path integral formulation of a quantum theory is expected 
to be mathematically ill defined. Other conceptual problems, such as the quantum 
effects at singularities, the problem of the vanishing of the cosmological constant, the 
meaning of measurements, are also far from being properly understood. Apart from 
the problem of the construction of a realistic model which could make contact with the 
low energy standard model and with cosmological observations, superstring theories 
are believed to solve the first problem, and to provide a consistent unitary quantum 
theory which contains gravity, at least in string perturbation theory. However, non- 
perturbative issues, such as the very existence of a vacuum state, the emergence of a 
realistic cosmology, the fate of black holes, etc., are still very open in these theories. 

In the lattice formulations of quantum gravity, one tries to construct a discretization 
procedure which makes the functional integration over metrics (or over the other 
degrees of freedom which describe classical gravity) meaningful, by truncating this 
infinite dimensional integral into a sum over a finite number of variables. Then 
one should look for critical values of the parameters (the coupling constants) of 
the discretized theory where a scaling behavior occurs, so that the details of the 
discretization become irrelevant and a continuum theory can be constructed. This 
lattice approach has been much successful for studying ordinary relativistic quantum 
field theories, such as non-abelian gauge theories, and in some cases to construct them 
rigorously. When trying to apply these ideas to General Relativity one should keep 
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in mind that one is making two basic assumptions: 

> One hopes that a quantum theory can be obtained from a path integral formulation 
based on the Lagrangian formulation of General Relativity, namely from the 
Einstein-Hilbert action, which is a functional of the metric tensor g, taken as the 
dynamical variable. There are some lattice formulations based on the first order 
Palatini formalism of GR, where the dynamical variables are the metric and the 
connection (or more adequately the local frame and the spin connection). I do 
not feel qualified enough to discuss these aspects here. I refer in particular to the 
lectures by A. Ashtekar for a thorough presentation of the Hamiltonian formulation 
of General Relativity and its canonical quantization. 

> One assumes that it is possible to recover Lorentzian gravity (where the metric 

has signature ( \- -\ — h)) from Euclidean gravity (where the metric has signature 

(+ + H — h)), by some kind of Wick rotation (this is the essence of the proposal by 
Hartle and Hawking (1983) for the wave function of the universe). 

Therefore, we are led to consider the functional integral over Euclidean metrics 



(D is the dimension of space-time, Gn Newton's constant and A the cosmological 
constant. I have set h = c = 1). 

One should also keep in mind that some very important problems have to be solved, 
and are not cured automatically by a lattice discretization. 

> Reparametrization invariance is the basic symmetry of GR. What does this 
symmetry become in a discretized theory, where space-time is not continuous? If it 
is broken by the discretization, how to ensure that it is restored (at least at large 
distances) in a continuum limit? 

> Are topology changing quantum fluctuations important, and how to deal with such 
fluctuations in the functional integral? 

> What about the unboundness of the Einstein action? 

> There is an important difference between a discretization procedure, which consists 
in making the number of dynamical variables finite or at least discrete, and a 
regularization procedure, which consists in introducing a fundamental energy scale 
(a regulator) above which quantum fluctuations are suppressed, in order to make 
the theory finite. Some regularizations (for instance higher derivatives theories) are 
not discretizations, on the other hand some discretizations cannot be considered as 
true regularizations, if they do not incorporate a non-zero minimal short distance 
cut-off. 




(1.1.1) 



with the Euclidean Einstein-Hilbert action 




(1.1.2) 



2. Piecewise flat manifolds, Regge calculus and dynamical triangulations 
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2.1. Simplicial manifolds and piecewise linear spaces 

Let me start with some geometry. I refer to Spivak (1970), Choquet-Bruhat et al. 
(1982) and to Dubrovin et al. (1984) for background in differential and Riemannian 
geometry. A simple and natural way to discretize space is to use lattices which are 
build out of simplices. I recall that a (oriented) (/-dimensional simplex is made out of 
d + 1 points pi, . . . Pd+ii the vertices, modulo (positive signature) permutations. For 
instance, a 0-simplex is a point, a 1-simplex a pair of points, i.e. a line or an edge, 
a 2-simplex a triangle, or a face, a 3-simplex a tetrahedron, etc. . . A subsimplex of a 
given simplex is simply a subset of points of this simplex. 




Fig. 1. 0-3 dimensional simplices, simplicial complex 



Out of a collection of simplices one can construct a simplicial complex, by gluing 
simplices along some of their subsimplices. Simplicial complexes may seem quite 
trivial objects. However one can already define exterior calculus on them. A p-form / 
is simply a function which associates to each p-simplex a some number /(c) (or more 
generally some element of an Abelian group), and such that, if —a is the orientation- 
reversed simplex, f(—a) = —f(a). One can define exterior product of forms, and an 
exterior derivation operation d, which maps p-forms into (p + l)-forms, by 

d/((p + l)-simplex) = /(p-simplex) (2.1.1) 

p — subsimplices 
of the simplex 

with the correct orientation. For instance 

df(pi,P2) = ffa) ~ f(Pi) 
df(pi,P2,£»3) = /(pi,P2) + f{P2,Pz) + /(P3,Pi) • (2.1.2) 

d satisfies the usual rule d 2 = 0. Therefore one can define cohomology theory on a 
simplicial complex and apply Stokes theorem 

/ df = f f . (2.1.3) 

JD JdD 

From a physicist point of view, this means that one can define scalar fields (0-forms), 
gauge fields (1-forms), antisymmetric tensors, etc. . . on any simplicial complex (see for 
instance Sorkin (1975)). However, a simplicial complex does not look like a manifold 
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with definite dimension, since it may be constructed from simplices with different 
dimensions. 

A (/-dimensional simplicial manifold (SM) is a simplicial complex such that the 
neighborhood of each point, i.e. the set of simplices which contain this point as a 
vertex, is homeomorphic to the (/-dimensional ball We shall often denote SM as 
triangulations, even when they are not two dimensional. 




Fig. 2. One and two dimensional simplicial manifolds 

A (/-dimensional SM can be obtained by successively gluing pairs of (/-simplices 
along some of their (d— l)-faces, until one gets a complex with no boundary. However, 
when performing such an operation, one does not obtain in general a SM, but only 
a so-called pseudomanifold (PM). Indeed, the complex constructed in that way are 
such that some vertices have neighborhood which are not homeomorphic to An 
example of 3-dimensional PM is depicted on Fig. 3. 




Fig. 3. An example of three dimensional pseudomanifold 

It consists of 4 tetraedra, 8 triangles, 5 edges and 3 vertices, and is obtained 
from the "double pyramid" by identifying the opposite faces of each pyramid, i.e. 
UNW <-> UES, UWS <-> UNE, DNW <-> DES and DWS <-> DNE } so that the 
"base" NESW becomes a 2-torus. The boundary of the neighborhood of each summit 
(U or D) is made of the two base triangles, and is homeomorphic to the 2-torus T2, not 
to the 2-sphere S2 = dBs ! However, all one or two dimensional complexes obtained 
by gluing lines at end-points or triangles along edges are indeed SM. This fact makes 
discretized models for two dimensional gravity much easier to handle. 
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The simplicial manifold structure corresponds to a discretization of ordinary 
manifolds, but there are no notions of distance between points, nor of volume element, 
which would allow to define angles or scalar products between vector fields, as well as 
to integrate functions over the manifold. One can embody a SM with a Riemannian 
metric by specifying that: 

> The metric is flat inside each d-simplex (i.e. the curvature vanishes). 

> The metric is continuous when one crosses the faces ((d — l)-simplices). 

> Each face of a simplex is a linear flat subspace of this simplex (i.e. its extrinsic 
curvature also vanishes). 

One thus obtains a piecewise linear space (PL). The metric on a PL space is entirely 
specified by the distance elements between nearest neighbor vertices, i.e. by the 
squared length of every edge of the corresponding SM. For instance, if we 

consider a triangle with vertices (1,2,3) and with edges length /12, I23, hi } hi the 
coordinate system where the metric inside the triangle is constant and the coordinates 
of the vertices are respectively (0,0), (0, 1) and (1,0), the metric is 

9 = \ 2 2^ 2 2 (^12 + ^31 — ^23)/ 

\(l 12 + ^31 — ^23)/^ ^31 

The number of edges of a simplex, d(d + l)/2, equals the number of independent 
components of the metric tensor. The constraint that the metric must be Riemannian 
means that the metric tensor is a positive definite matrix. It implies "generalized 
triangle inequalities" between the More generally if one wants a metric with 

signature (p,q), some I 2 - may be negative so that g^ inside each simplex has p 
positive and q negative eigenvalues. 

Finally let me recall that, while the metric is continuous when crossing the faces, it 
is discontinuous across the lower dimensional simplices. In particular, the curvature 
is concentrated along the (d — 2)-dimensional "hinges" (Regge (1961), Misner et al. 
(1973)). 

2.2. The dual complex and volume elements 

On a PL space we can construct the so-called dual complex, which is a natural 
generalization of the Voronoi construction of the dual lattice for a random lattice in 
Euclidean space (see in particular Christ et al. (1982) and Friedberg et al. (1984)). An 
intuitive (but not fully rigorous) definition is as follows: To any given p-dimensional 
subsimplex a we associate the flat (d — p)-dimensional dual polytope a* that consists 
of all the points which are at equal distances of all the vertices of cr, but closer to 
these vertices than to the other vertices of the complex (in fact this definition has 
to be slightly extended to take into account "obtuse-angled" simplices, which give 
non-convex polytopes). 

Examples in two dimensions are depicted on Fig. 4. The dual of a triangle is a 
point, the dual of an edge is a segment, and the dual of a vertex is a polygon. The dual 
polytopes of the simplicial complex C form a polytopial (but not simplicial) complex, 
C* . Let V{(t) be the p-dimensional volume (in the PL metric) of the simplex cr, 



(2.1.4) 



T93/028 SIMPLICIAL QUANTUM GRAVITY AND RANDOM LATTICES 



7 




polygon dual of a vertex 





Fig. 4. simplicial manifold (filled) and the dual complex (dashed) 

and V(<7*) the algebraic (d — p)-dimensional volume of the dual polytope a* . These 
volumes are function of the metric, i.e. of the edges lengths. The (/-dimensional 
volume element associated both to a and to its dual a*) is simply 



A(a) = A(a*) = ^^V(a)V(a*) 



(2.2.1) 



It is the volume of the set of points closer to a than to the others simplices of C. 
These volume elements allow to define a global inner-product over p-forms (functions 
over p-simplices) by 



[a, b] 



p\(d — p)\ 



p — dimensional 
simplices a 



V(a) 



a(a)b(a) 



(2.2.2) 



and a * operation (Poincare duality) which maps p-forms a on C onto (d — p)-forms 
*a on the dual complex C* (i.e. functions over (d — p)-dimensional polytopes) by 



(2.2.3) 
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This allows to define in full generality actions for fields on a PL space. For instance, 
the massless free field action for a scalar field $ (a 0-form) can be written as [d$, d$] 
and becomes 

fdx^g^d^d^ - ]T . (2.2.4) 

J edges V'J ) 

It is important however to notice that the * operation is not a duality, like in the 
smooth case, since the dual complex C* differs from C. 

2.3. Spin connection, spinors 

In fact it is more natural and convenient to define matter fields on the dual complex 
than on the initial simplicial complex (see in particular Ren (1988)). This is needed 
for spinors fields, which require the introduction of local frames e a = (e^) which live 
on the tangent space, since there are no well-defined tangent and cotangent planes at 
the vertices of C, and more generally at any point of the (d — 2)-dimensional hinges. 
One can associate an orthonormal local frame e a {i) (a = l,...d) to each vertex i 
dual to a (/-simplex, such that 

e a {i)e b {i) = eZ(iy b (i)g^(i) = S ah . (2.3.1) 

One can parallel transport these frames across the (d— l)-dimensional faces by simply 
translating them along the dual links. This allows to define the spin connection u 
as the 2-form which to a dual link (ij) associates the rotation matrix io a i l {ij) = 
e a {i)^b{j) Then one can define covariant derivatives for spinor fields ty a {i) living on 
the sites of the dual complex, couple spinors to gauge fields living on dual links {ij), 
etc. . . 

2.4- Curvature and the Regge action (classical Regge calculus) 

The curvature can be defined easily by performing parallel transport along a loop 
£ of the dual lattice which encircles a given {d — 2)-dimensional hinge h. Taking 
the ordered product of the rotations u{ij) we obtain a rotation matrix lo{£) which 
corresponds only to a rotation in the plane orthogonal to the hinge h, and thus is 
characterized only by the deficit angle Sh, which can be obtained as 

S h = 2tt - Y, 6 ^ , (2-4.1) 

where the sum runs over all the cZ-simplices meeting along h, and 8h,a is the dihedral 
angle of a along h. An explicit expression for 8h,a involves the volumes of cr, of the 
two {d — l)-faces of a which share h, and of h, that we denote respectively Vd, Vd-i 



and Vj_ 1? and Vd-2- It reads 
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The Regge discretization of the Einstein action is (Regge (1961)) 

5e = J d d xy^R -+ S Kegge = ]T V(h)6 h , (2.4.3) 

hinges h 

and is, through the volumes of the simplices, faces and hinges, a function of the edge 
length variables Uj. When varying the Regge action with respects to one edge length 
/ e , one obtains the Regge form of the Einstein equations (see Misner et al. (1973) for 
details) 

6ftcotan(0 ft ) = . (2.4.4) 

hinges h 
opposite to e 

The cosmological constant term in the action is simply 

I d d x^ -+ ]T V(a) . (2.4.5) 

d — simplices 

This discretized formulation of General Relativity, usually denoted Regge Calculus, 
is both elegant and very convenient for the study of various aspects of classical General 
Relativity. Indeed, when one approximates a smooth Riemannian space by a sequence 
of PL spaces obtained, for instance, by taking finer and finer triangulations of this 
space, with edge lengths given by the geodesic distance between the vertices, the 
sequences of metrics, curvature distribution, etc. . ., can be shown to converge, in a 
suitable sense, toward the smooth continuum ones. I refer to Christ et al. (1984), 
Cheeger et al. (1982) and (1984) and to Feinberg et al. (1984) for details on 
this classical continuum limit. Therefore, Regge calculus is an adequate framework 
to discuss classical gravitation (construction of classical solutions, discussion of 
gauge fixing). Small deformations of the edges length allow to deal with quantum 
fluctuations around classical solutions and to consider some semiclassical issues of 
quantum gravity. I refer to the recent and very complete review by Williams and 
Tuckey (1992) for references on the foundations of Regge calculus, recent advances 
and applications. 

2.5. Topological invariants 

Before addressing the problems of quantization, let me discuss shortly the issue of 
topological invariants for PL spaces. It is easy to write the Euler characteristic x, 
which is a non-trivial topological invariant for even-dimensional spaces, as 



X = Y^( — l) p N p ; Ni = number of p-simplices . (2.5.1) 

p=0 
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However, it is in general difficult to write simply other topological invariants, even 
when they can be written in terms of differential forms. For instance the Hirzebruch 
signature r for 4-dimensional manifolds is given by 

T = 7^ I d A x^e^°R ab ^R ab pa , (2.5.2) 

y(:,7r J M 

but for a PL space this involves products of distributions with support on the 
hinges. There is no simple regularization of these products which respects the 
topological invariance of r. I refer to Cheeger et al. (1984) for a general and 
detailed discussion of the definition and properties of curvature and of the simplicial 
approximation of manifolds by PL manifolds. 

2.6. Simplicial Euclidean quantum gravity (quantum Regge calculus) 

Can one use such discretization procedures to define the functional integral (1.1.1) over 
Euclidean metrics and to construct a quantum theory of Euclidean gravity? We have 
seen that it is possible to discretize the action. It is a more subtle issue to discretize 
the measure T>m [g] over the so-called "superspace" of all Riemannian metrics over 
some manifold Ai . In fact, since the local symmetry of GR is reparametrization 
invariance, one has to gauge away the diffeomorphism group Diffx, and to sum over 
orbits of Diffx in the superspace, namely to define J T>m M/Diff a4 . In the continuum 
a possibility is to start from the metric in superspace (see for instance DeWitt (1984)) 
given by the distance element 

IIMI 2 = / d d x^Sg„ [ag»»g vt > + Wg<">}8g p(T , (2.6.1) 
J M 

which defines a measure over metrics which is local and invariant under the action 
of Diffx, then to gauge away Diffx by using one's favorite gauge. Formally, this 
measure can be written as a product over all points of Ai (DeWitt (1967), Fujikawa 
(1983)) 

r" (d-4)(d + l) -r-f 

v\g] = n n d ^ > ( 2 - 6 - 2 

X fl<V 

but other choices of measure have been proposed, for instance with a different power 
for \g\ (Misner (1957), Leutwyler (1964), Fradkin and Vilkoviski (1973)). 

In order to define a discretized integration prescription over metrics out of Regge 
calculus, a first approach is to choose a given simplicial manifold C, and to consider 
all the admissible metrics (i.e. all the possible PL spaces) that one can construct out 
of C by varying the edge lengths l e . The integration over the / e 's will be restricted by 
the triangular inequalities which imply that the metric is positive definite on all the 
simplices. Then it is expected that reparametrization invariance will be restored in a 
continuum limit because one is summing over all the different piecewise flat metrics 
on C which are (at least approximately) equivalent by diffeomorphisms. 
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This approach, that we shall call here the fixed triangulation (FT) approach, has 
been mostly developed by Hamber and Williams (see Hamber (1986) for the basic 
theory, and Hamber and Williams (1984), and (1986)). In my opinion it has some 
appealing aspects but also some drawbacks. 

In the FT approach, many choices of measures for the edge lengths are possible. 
For instance, one can use simply a factorized measure of the form 

A J 

II -fVde) 7 > (2-6-3) 

edges 6 

e 

where V(/*) is the (d— 1)- volume of the face dual to / e , and f3 and 7 are some powers. 
More elaborate choice are possible, which involve for instance the volumes of the d- 
simplices. Let me mention that the natural expression in terms of edge lengths of the 
measure derived from the supermetric (2.6.1) involves a very non-local determinant 
depending on all the edges of the lattice. Another problem is that there can be in 
general very singular simplices, such as the triangles depicted on Fig. 5, with very 
small or very long edges. One usually deals with that problem by putting upper and 
lower bounds on the length variables. 



Fig. 5. Singular triangles 

However, such configurations are unavoidable when on tries to obtain a metric with 
very high curvature out of the flat metric. This is an example of the occurrence of 
what is is my opinion the main problem in this approach. It is not clear whether one 
can explore by this process regions in the space of metrics where the metric is very 
different from the typical metric on the reference simplicial complex that we have 
chosen at the beginning. 

These problems are related to the issue of gauge fixing and of the meaning of 
diffeomorphism invariance on simplicial spaces. This is discussed in particular in 
Frohlich (1982), Rocek and Williams (1984) and in Lehto et al. (1986). 

2.7. Dynamical triangulations (DT) 

A possible solution to these problems is to extend the sum over metrics by summing 
also over all topologically equivalent simplicial manifolds (if one fixes the topology), 
or even by summing over all possible simplicial manifolds (if one wants to deal 
with arbitrary topologies). A drastic simplification is then to reduce the sum over 
admissible metrics on a given C to only one representative metric, for instance the 
metric such that all edge lengths are equal to some fundamental length a. This length 
a will correspond to the ultra-violet cut-off of the theory. This approach has been 
first proposed for two dimensional gravity (see Frohlich (1985), David (1985a) and 
Kazakov (1985)), and then extended to higher dimensions. For some previous related 
ideas, see Weingarten (1977) and (1982). In the DT approach, one replaces the infinite 
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sum over the infinite dimensional space of metrics by a discrete sum over simplicial 
manifolds, which is mathematically more manageable. 

/ V M [g]/Ditf M ]T . (2.20) 

J SM 

For any "regular enough" Riemannian space, one can construct a sequence of PL 
spaces by taking a regularly distributed set of points, constructing a lattice out of these 
points and then taking as edge length the geodesic distance between nearest neighbors. 
Thus we can expect that the discrete set of PL space that we are considering is 
"regularly" distributed in the space of all metrics (Romer and Zahringer (1986)), 
while FT are localized in a neighborhood of smooth metrics on the regular lattice 
that has been chosen. Moreover we have introduced a natural UV cut-off in this 
approach. 




(a) (b) 

Fig. 6. Fixed triangulations with varying edges lengths are good approximations of metrics 
close to the flat one (a), but not of some other metrics, even when the curvature is bounded 
by the short distance cut-off (b) 

On the other hand we have completely lost diffeomorphism invariance, since one 
cannot deform smoothly (even approximately) a simplicial lattice into another one. 
One has to check if diffeomorphism invariance is recovered in the continuum limit, 
if it exists. Finally, it is not clear how to recover a classical continuum limit in this 
approach. 

Another problem, both conceptual and practical, for instance in numerical sim- 
ulation, is how to construct topologically equivalent triangulations. In fact such a 
notion of topological equivalence is a non-trivial one, and corresponds to (at least) 
two different concepts of equivalence. Two simplicial manifolds are homeomorphi- 
cally equivalent if they are equivalent by a (in general non-smooth) diffeomorphism. 
They are combinatorially equivalent if the two triangulations can be subdivised into 
the same finer triangulation. For two and three dimensions these two concepts are 
equivalent, but not in dimension d > 4 ! However for smooth triangulations they are 
always equivalent. I refer to Gross and Varsted (1992) for more details. 

It is better to consider only combinatorially equivalent triangulations, because they 
can be constructed by the so-called Alexander moves (Alexander (1930)), that I 
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describe now. Given a triangulation, take an edge e = (ij), then add a point k 
on e and subdivise the simplices which have e as edge into two simplices with edges 
(ik) and (kj) respectively. It is proven that one can always pass from one triangulation 
to another combinatorially equivalent triangulation by a finite series of such moves, 
with their inverses. 




2 dim 3 dim 



Fig. 7. Alexander moves in 2 and 3 dimensions. 

In practice another set of moves is often used, called (p, q) moves in Gross and 
Varsted (1992). Let C be a (/-dimensional triangulation. The boundary (i.e. the set 
of faces) of the ((/ + 1) dimensional simplex is a (/-dimensional complex, , which is 
a triangulation of the (/-dimensional sphere Sd- Take a subset S of p (/-dimensional 
simplices of C which can be mapped into a part of this complex and replace it by 
its complement S (it consists of q = d -\- 2 — p simplices) in E^. 




(1,4) (2,3) 

Fig. 8. (p,q) moves in 2 and 3 dimensions. 

It is known that such moves are equivalent to Alexander moves in 2, 3 and 4 
dimensions. Recently a general argument in any dimensions has been given by Pachner 
(Pachner (1991)). 
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3. Two dimensional gravity, dynamical triangulations and matrix models 

In this second lecture I shall discuss various formulations of Euclidean two dimensional 
gravity. Besides their interest as toy models for quantum gravity in higher dimensions, 
some models of two dimensional gravity, when coupled to two dimensional field 
theories, are in fact equivalent to models of string theories moving in some special 
backgrounds. This provides a strong motivation for their study, and explains some of 
their fascinating mathematical properties. I refer to the lectures by L. Alvarez-Gaume 
for more details. 

3.1. Continuum formulation of 2D gravity 

Classically, Einstein theory is meaningless in two dimensions, since the Einstein 
equations R^ v — ^g^R = are automatically satisfied. Indeed, the Einstein-Hilbert 
action is a topological invariant, since for closed orientable 2d manifolds we have the 
Gauss-Bonnet theorem 

/ d 2 xR = 4vrx = 8vr(l-/i) , (3.1.1) 

J M 

where \ is the Euler characteristic (2.5.1) and h the genus of Ai (number of handles). 
Nevertheless, some arguments indicate that a non-trivial quantum theory does exist. 

> Early renormalization group calculations by Weinberg (1979), following Christensen 
and Duff (1978) and Gastmans et al. (1978), for the theory in 2 + e dimensions 
(this makes sense by dimensional regularization, at least in perturbation theory) 
show that the scale invariance (and in fact the whole conformal invariance) of two 
dimensional gravity is spoiled by quantum anomalies. But if the quantum effective 
action is not scale invariant it cannot be a topological invariant ! 

> Polyakov (1981) stressed the crucial importance of conformal anomalies for string 
theories and 2d gravity. Starting from the functional integral (1.1.1) with the 
classical action (1.1.2) (with a cosmological constant), he showed that in the 
conformal gauge, i.e. in isothermal coordinates where 



9„u(x) = e**(*) , (3.1.2) 

the quantum action for the so-called "Liouville field" $^(x), obtained by integrating 
out the other matter fields and by taking into account the gauge fixing Fadeev- 
Popov determinant, can be calculated explicitly from the conformal anomaly and 
is still local in $£. It is the celebrated Liouville action 



26 



48vr 



d 2 ; 



\{d<S> L f + A R e** 



(3.1.3) 



where the number c is the conformal charge of the matter sector (it measures the 
amount of violation of scale invariance induced by the matter fields). Ar is the 
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renormalized cosmological constant, obtained from the bare one, A, by some UV 
divergent vacuum energy subtraction. 
> Finally, as discussed by Teitelboim (1983), the following classical equation for the 
metric (constant curvature equation) 



R + A R = (3.1.4) 

makes sense in two dimensions, and can be used as a starting point for quantization. 
It is nothing but the classical Liouville equation, 



-A$ L + A R e* £ = (3.1.5) 

(hence the name Liouville field for $^) obtained by the action principle from the 
Liouville action, rewritten in arbitrary coordinates. 

The conformal gauge is a very convenient framework to quantize 2d gravity, 
although other choices, such as the light-cone gauge, have some advantages (see 
Polyakov (1987) and Knizhnik et al. (1988)), in particular to bring forward the 
underlying current algebra of the model. The classical Liouville equation is integrable, 
and the program of the quantization of the Liouville theory has been initiated by 
Curtright and Thorn (1982) and by Gervais and Neveu (1982), and pursued actively 
by these last two authors and their collaborators (Gervais and Neveu (1983-85), 
Gervais (1991)). Many interesting results can already be drawn from the massless 
Liouville theory at Ar = 0, which is a free field theory, but is modified by boundary 
terms that I have not written in (3.1.3), which correspond to adding a "background 
charge" at infinity. I refer to the lectures by G. Zuckermann and L. Alvarez-Gaume 
for more details. 

An important result concerns the scaling properties for conformal field theories 
coupled to 2d gravity, and that of the partition function itself (Knizhnik et al. 
(1988), David (1988), Distler and Kawai (1989)). Let us consider the functional 
integral (1.1.1) for 2d gravity, with the action (1.1.2), where the metric g^ is now 
coupled to some scale invariant field theories, described by some fields that we 
call matter, with an action S'mf^I/, g^v}- The generating functional for field operators 
Oi(x) integrated over the whole space, Oi = J d 2 x ^J~g Oi(x), is obtained by adding 
source terms of the form JiOi to the action, 



Z[G N ,A,Ji] = 




(3.1.6) 



The cosmological constant A is proportional to the source Jo for the operator unity 
Oo = J <i 2 x y/gl, and we recast it in the J's. From (3.1.1) Z can be rewritten as a 
sum over the contribution of manifolds with fixed topology (given by the genus h) as 

oo l — h 

Z[G n ,Ji] = ( e ^) z h[Ji] • (3- 1 - 7 ) 

h=0 
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By taking derivatives w.r.t. the J's we define correlation functions (i.e. vacuum 
expectation values) for local operators averaged on a 2d space with fixed genus h } but 
with fluctuating metric, by 

(3.1.8) 

Ji = 



(0 1 ...O p ) h 



d 



d 



dJ\ ' ' dJ f 



Z h [Ji] 




Fig. 9. Correlation function for (microscopic) local operators 



If the metric g^ did not fluctuate, under a global dilation g^ 
correlators would scale with the total area A = J d 2 yjg 



as 



Xguv these 



{Ox... O p ) oc A-i 



(3.1.9) 



where A® is the scaling dimension (or conformal weight) of the local operator Oi(x). 
When the metric is quantized, the area A is no more fixed, but one can show that its 
average value (A) scales with the renormalized cosmological constant (see eq. (3.1.3)) 
as (A) oc 1/Ar. Then the scaling law given by eq. (3.1.9) is modified into 



£(i+4,-)+(2-7.)(i-ft) 



{O x ...O v ) h oc A 



R 



(3.1.10) 



where A{ is now the scaling dimension of the operator Oi "dressed" by gravity. The 
additional term in the exponent, proportional to the Euler characteristic \ = 2(1 — h), 
can be viewed as a topology dependent divergent part of the vacuum energy caused 
by the quantum fluctuations of the metric. The exponent 7^, often denoted the string 
exponent, depends only on the matter central charge c through the equation 



1 -7. + 



1 



13 



(3.1.11) 
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(valid for unitary theories) and is given by 

7a = -L ( c - 1 - ^(25 - c)(l - c)) . (3.1.12) 
The dressed scaling dimensions Ai are related to the bare dimensions A\ by 

Ai - A ^ l ~ A ^ = A° . (3.1.13) 

1 - 7* 

These scaling laws are often referred to as KPZ scaling, and can be derived easily 
from the Liouville theory. In the conformal gauge, one can show that the renormalized 
Liouville action for that has to be used in the functional integral given by eq. (3.1.6) 
is (taking into account quantum corrections to eq. (3.1.3)) 



S[*l] = j d 2 xyf§ (^gO'd^Ld^L + £$ L ) , (3.1.14) 

where g^ is a classical background metric and R the associated background curvature. 
The local operators Oi are "dressed" by gravity and become 

d = f d 2 x^/jjOi(x) e Qi * £ W . (3.1.15) 



The consistency requirement that the physics should not depend on the background 
classical metric g fixes the Liouville coupling constant 25 — c in eq. (3.1.14) and the 
renormalization factors a; to be 



1 

ai = — 
12 



(25 - c) - ^(25-cXl -c + 2AAI) 



(3.1.16) 



The dressed scaling dimensions Ai are given by the following formula, equivalent to 
eq. (3.1.13) 

^ = — , (3.1.17) 

a 

where «o is the renormalization factor for the operator with lowest dimension A 
coupled to the area, which is generically the dimension operator Oq = 1. 

It is much more difficult to compute explicitly the full correlation functions. Explicit 
results are known for genus h = and 1 (Gupta et al. (1990), Goulian and Li (1991), 
Di Francesco and Kutasov (1991), Dotsenko (1991). In addition to the correlations 
functions for local operators (sometimes denoted microscopic operators, one can also 
consider correlation functions for macroscopic operators W(£). Such an operator 
creates a loop with length £. In the presence of matter one must specify also the state 
of the matter along the loop. The corresponding macroscopic loop correlation 
functions (W{£\) . . . W(£ p ))h are obtained by summing over all metrics on a surface 
with h handles and p boundaries with fixed lengths £i } . . ., £ p . Such objects are much 
more difficult to compute in the framework of the Liouville theory. 
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Fig. 10. Correlation function for (macroscopic) loop operators 



3.2. Dynamical triangulations and continuum limit 



Following the ideas of dynamical triangulations developed in the first lecture (sect. 2), 
it is possible to define regularized models of two dimensional gravity by approximating 
the integration over Riemann surfaces by a sum over random triangulations (Frohlich 
(1985), David (1985), Kazakov (1985)). Restricting oneself at the moment to oriented 
triangulations, one defines the partition function of the model as 



1 



oriented 
triangulations 



'CO 



e 



a N + a 1 N 1 + a 2 N 2 



(3.2.1) 



where No, N\ and N2 are respectively the number of vertices, links and triangles of 
T, and s(T) is the order of the symmetry group of the triangulation T . ao, a\ and 
a>2 are free parameters. The symmetry factor s(T) is natural for the following reason. 
There are many different way to construct the same triangulation. We can start from 
sets of labeled vertices \y\ , . . . , ujv }, links {l\ , . . . , l^ 1 } and triangles {t\ , . . . , £jv 2 }, 
and glue them to get a labeled triangulation, and define the sum in eq. (3.2.1) as the 
sum over labeled triangulations, divided by iVoIiViIA^! which represents the number 
of labelings. Many labeled triangulations can be mapped one into another by a 
relabeling (i.e. a permutation) of their vertices, links and triangles, and therefore 
correspond to the same unlabeled triangulation. There are in general iVoIiViIA^! 
labeled triangulations for a unlabeled one, but for symmetric triangulations this 
number is smaller, and equals iVoIiViIA^! / ' s{T) , so that 

S NolN^l ^ 7CT) ' (3 ' 2 ' 2) 

labeled " " unlabeled ^ ^ 

triangulations triangulations 

For instance, the order of the triangulation made out of two triangles glued along their 
edges is s = 6, the order of the tetrahedron is s = 12. Since relabeling of elements 
of T can be thought as discrete diffeomorphisms, the factor 1/ s(T) can be viewed as 
what remains in the discrete case of the factor 1/Diffx. 
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From the topological relation \ = No — N\ + N 2 and the relation 2N\ = 3N 2 specific 
to triangulations, one sees that the discrete action in eq. (3.2.1) can be rewritten as 

a N + a 1 N 1 + a 2 N 2 = -\N 2 + 7X , 

A = - (j- + ^ + a 2 ) , 7 = a , (3.2.3) 

and that the action depends only on two independent parameters: A, that we identify 
with the cosmological constant A; and 7, that we identify with 1/Gn- Similarly we 
can define correlators for discrete loop operators Wk, (Wk 1 ■ ■ ■ Wk ), by the sum over 
triangulations with p boundaries, which are made out of k\ , . . . , k p links respectively. 

How to construct a continuum theory out of this discrete model? Let us first work 
with triangulations with fixed topology. The partition function Zh depends only on 
A. Then we introduce a fundamental length a (which corresponds to a short distance 
cut-off) that measures the size of the triangles. Up to inessential numerical factors, 
the dimensionful "physical quantities" are 

Area of the surface A = a 2 N 2 , 
Length of loop Li = aki , 
Cosmological constant A = A/a 2 . (3.2.4) 

To construct a continuum theory, we must find a critical point A c where the partition 
function is singular and the average number of triangles (N 2 ) diverges. Indeed, let us 
assume that at A c , Zh behaves as 

Z h (X) oc (A-A c )" a . (3.2.5) 

Then 

d a 

(N 2 )h = — "tt" In Zh = ^— + less singular terms . (3.2.6) 

If we define the renormalized cosmological constant A# as 

A R = a" 2 (A-A c ) , (3.2.7) 

we can define a continuum limit, where the short distance cut-off goes to zero, while 
physical quantities are kept finite 

a -> , A -> A c , A R fixed . (3.2.8) 

Indeed, in this limit, the average physical area (A), as well as its moments (A n ) } are 
finite 

(A) = a 2 (N 2 ) oc . (3.2.9) 

Ar 



T93/028 



F. DAVID 



20 



Similarly, one expects that the loop correlation functions, when expressed in terms of 
the physical lengths Li, are finite (up to possible wave function renormalizations). 

Of course at that point we have not said anything specific about the discrete model. 
In order to test these ideas one must 

> Solve analytically or numerically the model, and look for its critical point(s). 

> Check that the observables of the model are singular, and obey scaling relations 
which allow to define a continuum limit. In particular one should test universality, 
i.e. the irrelevance of the details of the regularization (one can use different lattices 
with various triangles, squares, add higher curvature terms to the action, etc. . .) in 
the continuum limit. 

> Finally, one should test the predictions of the continuum theory, in particular KPZ 
scaling, eq. (3.1.11-13), in order to check that the continuum theory that we have 
constructed is indeed two dimensional quantum gravity. In particular the exponent 
a in eq. (3.2.5) should be equal to (7^ — 2)(1 — h). 

3.3. The one-matrix model 

We have seen that to every triangulation T corresponds a dual simplex T*, which 
in the two dimensional case is made of polygons sharing trivalent vertices. Thus 
T* can be viewed as a Feynman diagram of a $ 3 -like theory. The basic idea of 
the matrix model formulation is to construct precisely a field theory whose Feynman 
amplitudes are the contributions of the triangulations in eq. (3.2.1). It is inspired 
from the ideas of 't Hooft (1974) about the large N c (number of colors) limit and the 
topological expansion for Quantum Chromo-Dynamics, and of subsequent works on 
planar diagrams and the planar limit (Brezin et al. (1978)). 

Let us consider the following matrix integral over N X N Hermitian matrices $ 
($ = $t) 

Z(g) = J d w2 $ e-™^*)) , (3.3.1) 

with the polynomial "potential" 

= -$ 2 - -$ 3 , (3.3.2) 
2 3 

and the flat measure 

d$ = C N JJ d$ 4J JJ dRe($ ij ) dlm^) , (3.3.3) 

i i<j 

normalized by the constant Cjy so that Z(0) = 1. This integral is not convergent, and 
at this stage we just consider it as a compact representation for its series expansion in 
powers of g, whose terms are well-defined. The term of order g k is given by the average 
([Tr($ 3 )] fc ) w.r.t. the Gaussian measure </$ exp(-f Tr($ 2 )), and can be expressed 
in terms of Feynman diagrams by applying Wick's theorem. The propagator can 
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be represented in terms of a double arrow (each arrow flowing from line to column 
indices) as 

($ij$ki) = jjSiiSjk J * == * • (3.3.4) 

There is no momentum flowing through the propagator because we are dealing with 
a simple integral, i.e. a "zero dimensional field theory". The cubic interaction term 
is 



gN$ij $ jk $*i 




(3.3.5) 



and the Feynman diagrams are trivalent "fat graphs" such as the ones depicted on 
Fig. 11. 




Fig. 11. a diagram and the dual triangulation 



The generating functional for the connected vacuum diagrams is F = ln(Z), and the 
contribution of a diagram Q is, counting the factors N coming from the propagators, 
the vertices, and the sums over indices which flow along each closed oriented loop, 

1 g n„ N n„-n p + n, ^ (3.3.6) 



where n„, n p and n\ are respectively the number of vertices, of propagators and 
of loops of the graph (?, while s(Q) is its symmetry factor. It is easy to see that 
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any connected Feynman diagram (fat graph) generated by F is the dual T* of some 
oriented triangulation T of the discrete model of subsec. 3.2. Since we have n v = ri2 } 
n p = ni, n\ = no and s(Q) = s(l~), we can identify the partition function Zdt of the 
DT model (eq. (3.2.1)) with the generating functional F for connected diagrams of 
the matrix model 

Z DT (A, 7 ) = ]n(Z(g,N)) ; g = e~\ N = e~'f . (3.3.7) 

Similarly, one can rewrite the correlators for the loop operators Wk as the connected 
correlation functions for the observables Tr($ fc ) 

{W kl ---W kp ) = (Tr($ fcl )---Tr($^)) conn . . (3.3.8) 

From eq. (3.3.6) a diagram with genus h contributes with a factor jV 2 ( fe_1 ). Therefore 
the topological expansion (eq. (3.1.7)) for two dimensional gravity corresponds to the 
1/N 2 expansion for the matrix model. The so-called planar limit, where only the 
lowest genus contribution is retained, is obtained by taking the limit N — > oo. 

The matrix models provide a surprisingly powerful method to study two dimen- 
sional gravity and string theories. Many matrix models are exactly solvable, and 
share very interesting mathematical properties with other fields of mathematics and 
physics (topology, integrable systems, etc. . .). They are discussed in details in the 
lectures by L. Alvarez-Gaume and E. Brezin. I shall discuss only a few simple points 
here. Some more elaborate discussion about matrix models and two dimensional 
gravity beyond the topological expansion and non-perturbative effects will be given 
in sect. 5. 

The most powerful method to discuss the matrix models is to use the invariance of 
the matrix integral under unitary transformations 

$ -> U$U^ ■ UU^ = 1 , (3.3.9) 

to diagonalize the matrix, and to rewrite the integral over the N 2 matrix elements 
as an integral over the N real eigenvalues of $. There is another method to 
solve the model in the planar limit, which is also of some interest. It relies on 
the Schwinger-Dyson equations for the matrix models (called the loop equations), 
and on the factorization properties of the observables in the large N limit (Wadia 
(1981), Migdal (1983)). If one performs the infinitesimal analytic change of variable 
$—)•$ + e/($) in the integral given by eq. (3.3.1), one obtains the identity 

iV(Tr[/($)y'($)]) = j ^M^(^)M^)) ■ (3-3.10) 

The l.h.s. comes from the variation in the potential V($), and the r.h.s. from the 
variation in the measure c/$. Rather than using the loop operators Wk = Tr($ fc ), it 
is more convenient (but mathematically equivalent) to consider the operator 



W(L) = Tr(e L *) 



(3.3.11) 
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(L is a continuous "pseudo-length"). Indeed, if one takes f(z) = e Lz , eq. (3.3.10) 
becomes the loop equation 



N V'(Jp-) (W(L)) = f L dL' (W(L')W(L — L')) . (3.3.12) 

Since the topological expansion of connected correlation functions of p loop operators 
involves surfaces with Euler characteristics \ ^ 2 — p, they behave as N 2 ~ p . It is a 
general feature of the limit N — > oo that non-connected correlations functions factorize 

(W(Li) ■ ■ ■ W(Lp)) ~ {W(L 1 ))---{W(L P )) . (3.3.13) 

Thus in the planar limit the loop equation (3.3.12) becomes a simple integro- 
differential equation for the one-loop correlator. The loop equations can be also 
used to compute multi-loops correlators, by adding to the potential V($) a source 
term J{L) for the loops 

r-OO 

— ► - / dL J(L) e L * . (3.3.14) 

Jo 

One obtains 

poo 

+ / dL" L" J(L")(W(L + L")) j . (3.3.15) 
Jo 

By taking derivatives w.r.t. J{L) one obtains multi-loops correlators. The first term 
corresponds to the splitting of the loop with "pseudo-length" L into two loops, while 
the second term corresponds to the fusion of the loop L with a loop L" . Thus the 
loop equations have a simple geometric interpretation in terms of evolution in loop 
space (see Fig. 12). They correspond to constraints for the partition function which 
form a Virasoro algebra. 

The solution of the matrix model in the planar limit (N — > oo), obtained either by 
the integration over eigenvalues or by the loop equation, shows that the discretized 
model of dynamical triangulations has indeed a critical point, where one can construct 
a continuum limit. This continuum limit corresponds to "pure gravity". In particular 
the exponent a obtained through the matrix model is equal to a = —5/2, which 
agrees with the KPZ scaling prediction for pure gravity (c = 0). 

3.4- Various matrix models 

Let me give some examples of matrix models which describe more complicated 
geometries, or gravity coupled to various kind of matter fields. 

General potential 
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V 







l-V 



V 




l-V 



Fig. 12. Loop equation for the Hermitian matrix model 



The one-matrix model can be solved for a general potential V of the form 

1 $2 _ H $ 3 _ 9± $ 4 _ H $ 5 _ 

2 3 4 5"' 



(3.4.1) 



This corresponds in the random lattice formulation of gravity to sum not only over 
triangulations, but over lattices made of triangles, squares, pentagons, etc. . . (Kazakov 
(1989)). For positive coefficients g n , one can show that the continuum limit is the 
same for all these models, and still describes pure gravity. 

When some of the coefficients g n are negative, one can obtain different continuum 
limits (Kazakov (1989)). It has been shown that they correspond to coupling gravity 
to non-unitary matter fields (Staudacher (1990)). For instance, the simplest non- 
trivial model corresponds to the CFT which describes the Lee- Yang singularity for 
the Ising model in imaginary magnetic field (i</> 3 model). 

Symmetric matrix model 

Another interesting model consists in replacing in eq. (3.3.1) Hermitian complex ma- 
trices by symmetric real matrices. In that case the sum over orientable triangulations 
is simply replaced by the sum over all triangulations, including non-orientable ones 
(Brezin et al. (1978)). This can be seen explicitly in the Feynman diagrams of the 
model, or in the loop equations. For instance, eq. (3.3.12) is replaced by 



nv ' { Jl ) {w(l)) = \[ dL ' mLr) W(L ~ L ' ]) 
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+ -L(W(L)) . (3.4.2) 

The first term, as previously, represents cutting a loop into two loops. The new term, 
which is of order iV _1 , is an orientation changing term representing the cutting of a 
loop bounding a Moebius strip (see Fig. 13). 



V 




i i-r i 

v 




i i-v 

Fig. 13. Loop equations for the symmetric matrix model 
Multi-matrix models 

One can generalize this discussion to multi-matrix models. The simplest example 
is the two-matrix model (Kazakov (1986), Boulatov and Kazakov (1987)). It 
corresponds to an Ising model living on the random triangulation, i.e. to Ising spins 
a (taking values ±1) on each triangle (each vertex of the dual Feynman diagram), 
with nearest neighbors coupling. One can implement this system by considering two 
matrices $_|_ and $_ and a potential of the form 

V($ + ,$_) = «($+ + $^) + 2/3$+$_ - g (e h <5>\ + e~ h $1) , (3.4.3) 

where h represents the external magnetic field, while a and f3 are related to the 
coupling J (in the interaction — J a a' between spins on nearest neighbors) by 
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In flat two dimensional space, the Ising model at its critical point is equivalent to a 
theory of free fermions, with central charge c = 1/2. 

The two-matrix model given by eq. (3.4.3) can be solved by the orthogonal 
polynomial methods (see the lectures by E. Brezin), as well as a more general class of 
multi-matrix models, which includes chains of matrices (Kostov (1989b). This class 
contains models of two dimensional gravity coupled to unitary rational conformal field 
theories, which have central charge c < 1. 

Another class of model which can be solved exactly (to a somewhat lesser extent) 
is provided by the model with one-matrix $ and n matrices \I/ a , a = l,n, with the 
following action (Duplantier and Kostov (1988) and (1990), Kostov (1988a)). 



Integrating over the \E''s, it reduces to a one-matrix model with a more complicated 
interaction, 



S[$] = Tr N -$ 2 - |$ 3 - - Tr N 2 [Ln [l<g)l - A(l®$ + $<g)l)]] (3.4.6) 



which describes a gas of self-avoiding loops (with line tension A and fugacity n) on 
random lattices. For fugacity —2 < n < 2 this system can be critical, and corresponds 
to the so-called 0(n) models, which include the self-avoiding walk (n = 0), the 
Ising model (n = 1), and the rotator model which describes the Kosterlitz-Thouless 
transition (n = 2). 

d- dimensional bosonic string 

Although one can obtain exact results for large classes of matrix models, many models 
have not yet been solved. This is the case of the (/-dimensional matrix model, where 
$(X) depends on the space coordinate X 



This model should describe a Polyakov bosonic string in (/-dimensional Euclidean 
space (Ambj0rn et al. (1985), David (1985b), Kazakov et al. (1985)). It can be 
solved in the large N limit only in the special cases d = (the ordinary one-matrix 
model), d = —2 (it can be mapped onto the 0(n = —2) model described above, see 
Kazakov et al. (1985) and David (1985c)), and d = 1, which is the continuous limit of 
an infinite chain of matrices. The solution of the d = 1 model is quite interesting, since 
it can be interpreted as a string theory in a D = 1 + 1 dimensional background. The 
additional coordinate is provided by the Liouville field of the effective string action 
(3.1.3). If one changes the kinetic term in eq. (3.4.7) so that the usual propagator 
becomes a Gaussian one 




(3.4.5) 





(3.4.7) 



1 




(3.4.8) 



K 2 + 1 
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(this is expected not to change the critical properties of the model), this amounts 
to changing the weight of each triangulation T in the partition function Z given by 
eq. (3.2.1) by a factor 



Wd{T) 



-d/2 



(3.4.9) 



where C-j is the No X No connection matrix of the triangulation T (it corresponds 
almost to the discretized scalar Laplacian operator A = d* d discussed in subsec. 2.2), 
and det' means the product of the No — 1 non-zero eigenvalues of the matrix —C. 
Moreover one has an explicit expression for the determinant in eq. (3.4.9) in terms of 
the number of spanning trees of T, i.e. of connected subgraphs of the triangulation 
which contain all the vertices of T, but have no internal loops 

— — det'( — C-j) = ft of spanning trees of T . (3.4.10) 

This weight corresponds to integrate over all embedding of the triangulation T in 
"physical" (/-dimensional space. One assigns a position vector X(v) to each vertex v, 
the discretized action for the embedding is 

S[X] = J2(X(v)-X(v')f , (3.4.11) 

edges 



while from the rules of subsec. 2.2, the measure over X becomes 

V[X] -> Y[ d d X(v)c(v) d / 2 , (3.4.12) 

vertices 

where c(v) is the coordination number of the vertex v (Exercise: show how this term 
comes from the rules of subsec. 2.2 ). One has finally 

/ V[X] e~ s ^ oc JJ c(v) d l 2 w d (T) . (3.4.13) 



3.5. Numerical studies of 2D gravity 

Besides the exact results that can be obtained through the matrix models, one must 
rely mostly on approximate numerical schemes to study the dynamical triangulation 
models. However rigorous estimates have also been obtained, and some will be 
discussed in due time. As discussed in subsec. 3.2, one is interested in the structure 
of the phase diagram and the possible critical points for the DT model defined 
through the partition functions of the form (3.2.1). The sum over triangulations 
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makes sense for fixed topology, and if not otherwise specified, we shall deal with 
planar triangulations (genus h = 0). The partition function is of the form 

z(g) = 7tT) w ^ gN2 > t 3 - 5 - 1 ) 

planar ^ ' 

triangulations 

where N2 is the area of the triangulation T and the weight w depends on the matter 
fields which live on the triangulation. The matrix model results show that the number 
of triangulations with fixed topology grow with the area as (constant) N2 . Moreover the 
constant does not depend on the topology (genus) of the triangulation, but depends 
on the details of the discretization (Rivasseau (1991a)). This implies that eq. (3.5.1) 
is a convergent series for \g\ small enough, which defines Z(g) as an analytic function 
around the origin. This function will have singularities in the complex plane. We 
are usually interested in the smallest singularity go on the positive real axis. This 
singularity will allow to get a continuum theory, along the lines of subsec. 3.2 (see Fig. 
14). I now present some of the numerical approaches to study the critical behavior. 

Series expansions 

Usually go is the closest singularity from the origin, and the critical behavior of Z(g) 
at go 

Z(g) ~ z c (go — g) 2 u + regular or less singular terms (3.5.2) 

is therefore related to the large order behavior of its series expansion 

00 

Z(g) = Y, z *9 k ; z k = z c k'''- 3 go- k {l + o(k- 1 )) (3.5.3) 

A;=0 




Fig. 14. Singularities of the partition function in the complex plane 
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One can compute the first terms Zk of the expansion by constructing all the 
triangulations with a fixed number of vertices (this can be done with the loop 
equations), and by estimating their weight factor w. By such methods one can 
typically compute the first 10 ~ 15 terms of the expansion (David (1985b), David 
et al. (1987), Ambj0rn et al. (1986b)). Then one obtains estimates for the value of 
the critical point go and the critical exponent j s . However in general this is made 
difficult because the subdominant terms in eq. (3.5.3) are not small. One reason for 
this is that there might be other singularities in the complex plane, beyond go } but 
close to it. Various recipes exist to deal with these problems. One can take adequate 
combinations of the Zk which depend weakly of these terms (ratio methods). One 
can also approximate Z(g), or rather some derivative of its logarithm ln(Z(^)), by 
rational fractions (d-log Pade methods). These kinds of analysis deal less easily with 
subleading singularities (such as (go — g) 2 1 , 7' > 7«) which occurs at the same 
critical point (confluent singularities). These singularities are expected to occur for 
DT describing some matter fields coupled to gravity, since they will correspond to 
perturbation by matter operators less relevant than the identity operator 1 coupled 
to the area via the cosmological constant. 

Monte-Carlo simulations 

Instead of enumerating exactly all the triangulations with a given number of triangles, 
one can explore the space of triangulations by Monte-Carlo simulations (Kazakov et 
al. (1986)). One constructs a stochastic process, which is defined by its transition 
probabilities: if at "time" t the system is in a configuration T, it has a probability 
P(T 1— > T) to move to configuration T at time t + 1. To reconstruct the sum over 
triangulations given by eq. (3.5.1), one must ensure the fact that in the large t limit, 
the probability to visit a configuration T is proportional to its weight w(T) (we have 
recast the factor g N2 in w) in the partition function. This will be the case if: 

> the process is ergodic, i.e. all configurations can be reached starting from an 
arbitrary initial configuration; 

> the probabilities P satisfy the detailed balance equations 

W(T)J2 p ( T ^ T ') = ^2W(T')P(T' h->T) • (3.5.4) 

T T 

We have discussed in subsect. 2.7 examples of moves (Alexander moves and (p, q) 
moves) which are ergodic in the space of (combinatorially) equivalent triangulations. 
In two dimensions the situation is simpler, since the "flip" (the (2, 2) move) can be 
shown to be ergodic in the subspace of triangulations with fixed number of triangles. 
Therefore one can perform M.C. simulations in the two following ensembles: 

> By using only the flip, one deals with the "canonical ensemble", where the number 
of triangles N2 (the volume of the manifold) is fixed. The limit N2 — > 00 corresponds 
simply to take the critical limit Ar — > 0. 

> By using the (1,3) and (3,1) moves (insertion and deletion of a vertex inside a 
triangle) one deals with the "grand-canonical ensemble", where the volume varies. 
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One has to fine tune the coupling constant g (the bare cosmological constant) to 

reach the critical point. 
When matter fields are present (for instance the Ising spins of the model described 
by eq. (3.4.3), or the d Gaussian fields X(v) = (X l (v); i = l,d) which describe the 
position of the vertices v of the triangulation in the external space for the model 
given by eq. (3.4.7-9)), the average over these fields can be performed independently 
by standard Monte-Carlo methods. 

The canonical ensemble does not allow in principle to estimate the susceptibility 
exponent j s . In practice, the use of ensembles interpolating between the canonical 
and the grand-canonical ensemble is quite convenient. For instance one can allow the 
volume of the system to vary between two bounds iV m i n and iV max , and extract the 
value of the critical coupling and that of 7^ from the study of the average distribution 
of N2 between these two bounds, either by direct fits, or by using a maximum 
likelihood principle. 

Up to now the most thoroughly studied model (both by series expansions and by 
M.C. simulations) is the (/-dimensional string model corresponding to the Gaussian 
weight of eq. (3.4.9). I refer to Kazakov et al. (1986), Boulatov et al. (1986), Ambj0rn 
et al. (1986), Billoire and David (1986), Jurkiewicz et al. (1986a-b), David et al. 
(1987), Ambj0rn et al. (1987), and to the more recent reviews by Ambj0rn (1991) 
and Kawai (1992) for details. 

Other models have also been studied recently, in particular systems of n Ising 
models coupled to a random triangulation (only the case re = 1 is exactly solvable) 
(Brezin and Hikami (1992)). 

Sampling methods using exact results 

In the special case c = 0, the weight factor of a triangulation is basically 1 (up to 
symmetry factor), and the number of open planar triangulations (with the topology 
of a disc) with fixed area (number of triangles) and fixed boundary length (number 
of edges) is known exactly (from the matrix model exact solution). This allowed 
Agishtein and Migdal (1991a) to construct samplings of triangulations via a biased 
random process, with the correct probability distribution, and to go to very large 
triangulations (up to 10 7 triangles, i.e. two orders of magnitude larger than with 
ordinary M.C. Methods). Of course this is not useful for evaluating the critical 
exponent 7^, which is exactly know, but this allows to study the "internal geometry" 
of 2d quantum gravity, and to measure quantities such as the "internal Hausdorff 
dimension" of space, which are not (yet ?) accessible through exact solutions. This 
method has been cleverly extended to the c = —2 case (Kawamoto et al. (1992)). I 
shall discuss briefly some issues about the internal geometry of 2d gravity in subsec. 
3.7. 

2d quantum Regge calculus 

For a recent review on numerical simulations of 2d gravity using quantum Regge 
calculus (fixed triangulation and varying edge lengths) discussed in subsec. 2.6, I 
refer to (Hamber (1991a) and to Gross and Hamber (1991)). Numerical simulations 
have been performed with or without a higher derivative term (the square of the deficit 
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angle) which amounts to add a aR 2 term in the action. Estimates for the 7^ exponent 
for lattices with the topology of the torus (genus h = 1) give 7^ ~ 2.025 ± 0.025, in 
agreement with KPZ scaling. For the genus case, the results of Gross and Hamber 
(1991), although less accurate, agree also with KPZ scaling. 

As discussed above, I think that the problems of the restoration of coordinate 
invariance, and how FT take into account large fluctuations of the metric away from 
the regular flat metric are not fully clarified. In particular the result 7^ = 2 is 
quite universal for genus h = 1 surfaces, since it does not depend on the matter 
central charge c, and is even valid when the metric does not fluctuate (c = —00). 
More numerical studies for genus h = lattices are in my opinion needed to really 
understand the model. 

Studies of the phase diagram of a model of FT in (/-dimensional function 
of the R 2 coupling a, of the form of the measure (see eq. (2.6.3)), and of the dimension 
<f, have also been performed. They indicate a transition between a "smooth phase" 
for positive a and a "rough phase" for negative a. 

The simulations of Ising model on 2d Regge lattice indicate no drastic modifications 
of the critical exponents of the Ising transition when compared to the Ising model on a 
regular flat lattice (Gross and Hamber (1991). This is at variance with the predictions 
of KPZ scaling, and the results of the Ising model on dynamical triangulations, and 
is an argument for the fact that these models do not belong to the same universality 
class than the DT models. 

3.6. The c = 1 barrier 

The numerical methods described previously give excellent agreement with the exact 
solutions for DT models coupled to statistical models corresponding to c < 1 matter. 
I shall mostly discuss here the numerical results for systems where no exact solutions 
are available. 

It is clear from the KPZ scaling formulae (eq. (3.1.12-13)) that the continuum 
formulation of 2d gravity makes sense in the weak coupling phase c < 1, but that 
something happens for c > 1, since the critical exponent 7^ becomes complex for 
1 < c < 25. In fact one expects also that the theory will be ill-defined for c > 25, 
since the kinetic term for the Liouville field becomes negative (the 26 — c in eq. (3.1.3) is 
renormalized into 25 — c by quantum corrections). Before discussing numerical studies 
of this regime, let me first explain what is the present theoretical understanding of 
this question. 

Surfaces versus branched polymers 

A general argument given by Durhuus, Frohlich and Jonsson (Durhuus et al. (1984)) 
suggests that, if the string exponent for spherical topologies, 7^, is positive, branched- 
polymer-like configurations should dominate the functional integral over 2d, metrics. 
In fact their argument is rigorous (up to very plausible scaling hypothesis) only for a 
specific class of random surface models made of plaquettes on a lattice in bulk space, 
but heuristically it sounds quite reasonable for the random triangulation models. 
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A simple model of abstract branched polymer can be constructed out of links which 
can join themselves at vertices in all possible ways to form a connected tree. To each 
vertex we assign a statistical weight of the form 

g ■ /(c) ; c = number of links attached at the vertex. (3.6.1) 

g is a "vertex fugacity", and /(c) some connectivity dependent factor. The 
coordination number c can go from 1 (end-point) to oo. The partition function of the 
model is simply 

Z(g) = 9 # ° tvertlcea II M • (3-6.2) 

trees vertices 

Z(g) becomes singular at a critical point g C} where the average number of vertices 
diverges. Generically, if /(c) decreases fast enough as c -» oo, this critical point 
describes branched polymers; the average coordination number stays finite and the 
critical exponent 7^ (such that Z(g) oc (g c — g) 2 u ) is that of branched polymers 

Is = \ ■ (3.6.3) 

The basic idea of Durhuus et al. (1984) is that if one decomposes planar surfaces into 
"blobs" connected by small "bottlenecks" (such configurations are always present in 
the integral over metrics), one can show that if 7^ > 0, the contribution of each 
blob to the sum over metrics grows as g — > g c in such a way that the total number 
of blobs and bottlenecks will diverge before the average size of each blob becomes 
large. It follows that, close to the critical point, i.e. for surfaces with large area, the 
configuration space is dominated by branched polymers made out of blobs (with a 
finite area), rather than by 2d surfaces. In fact for the dynamical triangulation model 
discussed in subsec. 3.4, it is shown rigorously that if 7^ > 0, then 7^ < 1/2, so that 
the true branched polymer exponent 1/2 represents an upper bound. 

A simpler, but heuristic argument (David (1992)), consists in comparing the area 
dependence of the partition function for one surface 

Z(A) oc A 73-3 (3.6.4) 

with that of a pinched surface of the same area, made by gluing two marked surfaces 
of area A' and A — A' at the marked point. Since the partition function of a marked 
surface should scale as 

Z m ^ ed (A) oc A ■ A ls ~ 3 = A ls ~ 2 , (3.6.5) 
the partition function for the pinched surface scales as 

Z pinc UA) oc / dA'A n °- 2 (A-A'y*- 2 ~ A 2 73-3 . (3.6.6) 
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Comparing eq. (3.6.4) and eq. (3.6.6) one sees that pinches are entropically favorable 
if 7 S > 0. For 7 S < (c < 1) such pinches are not favorable. However, the statistics 
of pinches which create small "baby universes" with area A' <C A leads to interesting 
questions, and can be studied via matrix models and the Liouville theory (see Jain et 
al. (1992)). 

Tachyons for the bosonic string 

From eq. (3.1.13) j s is negative for c < 1 and vanishes at c = 1. This can be related 
to the problem of the tachyon in bosonic string theories (see in particular Seiberg 
(1991)). The (/-dimensional non-critical bosonic string described by the continuous 
action (3.4.10) (with c = d) corresponds to a string in D = d + 1 dimensions. 
Indeed, if the Liouville field is identified with an additional coordinate X° by 

$i = yj ' 2 ^ d X° ', so that X = (X°,X), the Liouville+matter action becomes 

J_ f rr fl 

4tt 



f d 2 xyfg~ Q^^Xd^X + RD(X°) 



D(X°) = f\^X° , (3.6.7) 

where D(X°) is a linear background dilaton field. The operator e lPX dressed by 
gravity (see eq. (3.1.15)) can be identified with a tachyon vertex operator (i.e. an 
operator which creates a spin state of the string) in D = d + 1 dimension with 
momentum P = (Po, P) by 

e iPx+ ap <s> L = e D(x°) e iPx+p x° _ (3.6.8) 



The component Po of the momentum P is 

Po = ^M2+(P)2 • M 2 = l_£ . (3. 6 .9) 

For d < 1, Po is real. This suggests to consider X° as a time-like coordinate, Po as the 
associated energy-like momentum component, and M as the mass of the spin ground 
state of the string. For d < 1 the ground state of the string is massive, for d = 1 it is 
massless, while for d > 1 it is really tachyonic. Thus the c = 1 barrier is nothing but 
the D > 2 tachyonic barrier of the critical bosonic string. Seiberg (Seiberg (1991)) 
has related the existence of tachyons for d > 1 to the so-called "macroscopic states": 
for d > 1 the insertion of the puncture operator 1 (i.e. the tachyon vertex operator at 
zero momentum P = 0) does not create a small, local, disturbance of the metric, but 
a large, macroscopic disturbance, because the distances do not become small when 
one approaches the insertion point. A related problem is that for c < 1 there is only 
a finite number of states for 2d gravity coupled to matter (this can be viewed as the 
reason why the c < 1 models are solvable), while for c > 1 there is an infinite number 
of states. 
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Spikes/wormholes picture 

Finally another point of view (but somehow related to the previous one), consists 
in looking at the importance of singular configurations with conical curvature 
singularities (spikes) in the functional integral over metrics (Forster (1987), Cates 
(1988)). In the Liouville formulation of 2d gravity such defects appear as point like 
charges, with a charge q proportional to the defect angle AO at the conical singularity. 
A simple argument, similar to the well-known Kosterlitz-Thouless argument for the 
2d Coulomb gas model, shows that curvature defects such that \A9\ is smaller then 
some critical angle A(c) are liberated. Explicit estimates indicate that as long as 
c < 1, A(c) < 2tt and this gas of curvature defects just contributes to the fluctuations 
of the metric, but for c > 1, A(c) > 2tt and singular spikes (infinitely long and thin 
cylinders) are liberated (see Cates (1988), Krzywicki (1990) for details). 

Numerical studies of the c > 1 phase 

I now briefly summarize the numerical results obtained so far. As already mentioned, 
numerical studies agree quite well with exact solutions which corresponds to — 2 < 
c < 1, i.e. — 1 < 7 S < 0. The "strong coupling regime" c > 1 is much more delicate 
to explore. Here I shall describe mainly the results of Ambj0rn et al. (1986a, b), 
Jurkiewicz et al. (1986a, b) and of David et al. (1987) on the (/-dimensional DT model 
described in subsec. 3.4, which corresponds to a non-critical string in (/-dimensional 
Euclidean space. To probe the effects of higher order curvature terms in the action, 
as well as the importance of the form of the measure given by eq. (3.4.11), one can 
add a term of the form FJ cf , where Ci is the coordination number of the vertex z, 

vertices 

so that the weight of a triangulation T becomes 

VdAT) = il ^ ( # t °ee 3 3P onT S )"" /2 ■ (3-6.10) 

vertices 

Expanding C{ around its average value (ci) = 6, and using Euler theorem, one sees 
that 

Y[cf ~ e - ^ T ti ( Si r 1 Y+- , (3.6.11) 

i 

so that a > will decrease the average squared curvature R 2 and flatten the 
triangulation (at least at short scales), while a < will increase R 2 and make the 
triangulation more "crumply" . 

This is indeed what is observed, both in MC simulations and series analysis. For 
c < 1 and a > one is in a "Liouville phase" where KPZ scaling holds and where 
7 S is negative and independent of a. For a < and large enough one enters a new 
"collapsed phase" where 7^ becomes very negative and a dependent. The existence 
of two different regimes is also observed for c > 1. For a > and large, one is in 
a "branched phase" with 7^ ~ .5, in qualitative agreement with branched polymer 
models. For a < and large, one is still in the collapsed phase. 
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Fig. 15. Qualitative picture of the phase diagram for the (/-dimensional DT model 

The nature of the branched and collapsed phases may be understood qualitatively 
when c is large by analytical estimates (Ambj0rn et al. (1986a, b)). One can show 
that in this limit, triangulations obtained by gluing together tetraedra along edges 
dominate the partition function, since they minimize det'( — C-j) in eq. (3.4.9). Such 
configurations correspond to branched polymers. Depending on whether a ^> —d/2 } 
or«< —d/2, two very different kinds of BP on maximize the weight Wd, a - 





Fig. 16. Gluing two tetraedra 



If a <C —d/2 } it has been shown that very singular configurations, where all the 
tetraedra are glued along a single link, are dominant. This explains the collapsed 
phase. The R 2 coupling is so large and negative that the metric collapses into a few 
very singular regions where the curvature is very large and negative, which correspond 
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to the points with coordination number of the same order of magnitude than the total 
number of points of the triangulation. The susceptibility exponent is found to vary 
linearly with a as 7, ~ 2a-\-d/2. Such a phase has no continuum limit interpretation, 
since all points are at a finite distance from each others, and the Hausdorff dimension 
of such a space is infinite. 

On the other hand, if a ^> —d/2 (this includes the "physical situation" a = 0), it 
is ordinary BP, with no vertices with large coordination number, which dominates, 
and 7 S = 1/2. This corresponds to the "branched phase". 

Unfortunately, this picture is only valid for large <i, and checked numerically in a 
convincing way for typically d > 10. A qualitative understanding of the c > 1 phase 
for c close to 1 is still lacking. Numerical estimates for 7^ show a smooth increase from 
to .5 as c goes from 1 to 00, however finite size effects seem important, even for the 
largest triangulations used up to now in MC simulations. The same conclusion can 
be drawn from series analysis from the (/-dimensional bosonic string model. Recent 
series analysis of Ising models coupled to DT by Brezin and Hikami (1992) corroborate 
these conclusions. Therefore it is not clear yet whether for c > 1, 2d gravity is in a 
branched polymer phase, with 7^ = 1/2, as long as it is not in a singular, collapsed 
phase, or if there might exist some intermediate non-trivial phases with < 7^ < 1/2. 
Such phases would be extremely interesting to characterize. It would also be worth 
investigating the nature of the transition between the branched and collapsed phases. 

3.7. Intrinsic geometry of 2D gravity 

Let me discuss briefly numerical studies of the intrinsic geometry of 2d gravity. 
In addition to the global observables, obtained by averaging over 2d space local 
(microscopic) operators or macroscopic loop operators, one may consider correlations 
involving the geodesic distance between points. The simplest of such observables is 
the volume V{R) of the ball (in 2d, the disk) with radius R } defined as the set of 
points x at geodesic distance fj^g ^ X ^ X Q ) < R from the center x$. This allows to define 
the internal Hausdorff dimension dn as 

{V(R}) — > R dH . (3.7.1) 

R — >oo 

In a non flat 2d geometry the disk is connected but does not need to be simply 
connected, so that the number of connected components of its boundary, n, might be 
larger than 1 (see fig. 17). If this number scales with the radius R of the disc, this 
defines a new branching dimension 6b by 



{n(R}) — > R 6b 

R — >oo 



(3.7.2) 
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Fig. 17. A non-simply-connected disk 

For flat geometry one has obviously djj = 2 and 8b = 0. For DT these quantities 
cannot be calculated exactly from the matrix models. One has to use to very 
large triangulations (~ 10 6 triangles !) in numerical simulations to see some scaling 
behavior, and some qualitative estimates have only been obtained for the cases c = 
(Agishtein and Migdal (1991a)) and c = —2 (Kawamoto et al. (1992), by using the 
sampling methods described in subsec. (3.5). For c = —2 the Hausdorff dimension djj 
seems much larger than 2 (c/jj ~ 3.5), and the branching dimension is positive 8b > 1. 
For c = 0, da and 8b are even larger, but the numerical evidences for scaling are less 
convincing. This shows that, at the scales considered in numerical simulations, the 
2d geometry for DT is already very "spiky" in the phase c < 1. However, one sees 
that, from these numerical results, DT seem more branched than "ordinary" branched 
polymers, which are characterized by djj = 2 and 8b = 1 (see David (1992)), and 
more similar to infinitely branched trees with very few end-points, such as the Cayley 
tree. 

Let me mention also that analytical estimates of these quantities through Liouville 
theory is a difficult problem (see David (1992)). Indeed, KPZ scaling holds for field 
theories close to their critical points coupled to gravity, and therefore can be viewed as 
involving massless particles, or at least particles with masses m oc L~ x <C a -1 of the 
order of the IR cut-off (L ~ \J~A~ is the size of the manifold), and much smaller than 
the UV regulator (a is the UV cut-off). On the contrary, the geodesic distance can 
be expressed in terms of the propagation of very massive particles (£ _1 <C m < a -1 ) 
in the fluctuating metric. The propagation of such massive particles is very different 
from that of massless ones, and new physical phenomena, such as "trapping" by 
regions with high curvature, occur, which have to be taken into account properly. 
Interestingly, these problems are related to problems encountered in the statistical 
mechanics of disordered systems (spin glasses, propagation in random potentials). 

3.8. Liouville theory at c > 25 

There is a very important difference between the models of 2d gravity in the weak 
coupling phase c < 1 that we have discussed insofar, and Einstein theory of gravity in 
dimensions d > 2. For c < 1 the effective theory is the Liouville theory, whose action 
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is given by eqs. (3.1.3,14), with a positive coupling constant. This action is bounded 
from below, and for positive cosmological constant Ar (i.e. positive vacuum energy 
density), the stable ground state corresponds to anti-de-Sitter space (R < 0). 

Following Polchinski (1989), one can try to study the Liouville theory with a 
negative coupling constant, which corresponds formally to the case c > 25, as a toy 
model for the d > 2 Einstein theory with positive Newton's constant. The Liouville 
action, which takes a local form for conformally flat metric, can then be rewritten for 
general metrics in a bilocal form 



Because of the minus sign in the r.h.s. of eq. (3.8.1), this model admits de Sitter space 
R = Ar as classical solution for Ar > 0, but suffers from the conformal instability 
problem. 

One can study the problem of the stability of this action by the dynamical 
triangulation regularization. Since the configuration space is discrete (for fixed 
area), the action is in fact bounded from below. One might ask what are the real 
configurations which minimize this bilocal action, i.e. which maximize the Liouville 
action. Simple, but not rigorous arguments suggest that in that case also, branched 
polymer configurations minimize the action. Indeed, if one considers a genus metric 
with fixed area A, the Liouville action of the sphere scales as ln(A), while the action 
of a long cylinder with radius a (the UV cut-off) and length L oc A/a scales as 
—L/a. Moreover, one can check that the action for tree-like configurations made out 
of cylinders with width a is still proportional to —L/a } where L is the total length of 
the tree. Thus one sees that the action is still extensive, i.e. proportional to the area. 
Moreover, one can argue that such branched polymer configurations do extremize the 
Liouville action (David, unpublished). 




(3.8.1) 



C 




L 



Fig. 18. Configurations which maximize the Liouville action 



4. Euclidean quantum gravity at D > 2 



In this lecture I shall give a short review of recent results on lattice models for quantum 
gravity in dimensions D > 2. Since I have not worked myself on these models, my 
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review will be probably incomplete and will reflect some of my prejudices. I shall 
deal neither with the formulations of 3D gravity as Chern- Simon theories (Witten 
(1989)), nor with the 3D topological theories a la Ponzanno-Regge (Ponzanno and 
Regge (1968), Turaev and Viro (1992)), which allow to construct new topological 
invariants for 3d manifolds. Let me first discuss what one could expect from the 
lattice approaches of the problem of the quantization of the Einstein theory. 

4-1. What are we looking for? 

For D > 2, the Einstein action given by eq. (1.1.2) is non-renormalizable, since the 
Newton's constant is dimensionful, with Gn ~ [length] D ~ 2 . If one still wants to 
describe gravity by some local field theory, one possibility is to add higher derivative 
terms to the action. In D = 4, the curvature squared theory obtained by adding to 
the action the two independent dimension 4 operators 

d D Xy ^(aR 2 + f3 W 2 ) , (4.1.1) 

where W is the Weyl curvature tensor, is known to be renormalizable and asymp- 
totically free (Fradkin and Tseytlin (1982)). Thus it might define a Euclidean theory 
sensible at short distance, with dynamical generation of a mass scale (as in QCD), 
with General Relativity as a low energy effective theory below this scale. Moreover, 
the action is now bounded from below. However, there are still serious problems 
with this approach. In particular, such higher derivative theories are known to be 
non-unitary. Since the conformal instability problem is still present at large scales, it 
is not clear whether the effective low energy theory is really the Einstein theory at 
scales i > Lpi anck . 

Another approach has been advocated some time ago by S. Weinberg (Weinberg 
(1979)). From dimensional analysis, the effective coupling constant of Einstein theory 
is known to grow with the momentum scale p as G(p) oc p D ~ 2 . It is possible to study in 
perturbation theory the renormalization group flow of G(p) for space-time dimension 
D > 2 within a D = 2 + e expansion. Such an expansion makes sense in perturbation 
theory (see Weinberg (1979) and Kawai and Ninomiya (1990)). The interesting result 
is that the beta function, calculated at one loop, is of the form 

BC 

p— = /3(G) = eG - cG 2 + ... (4.1.2) 

dp 

with c some positive constant (at least when no matter fields are coupled to gravity). 
The standard interpretation of this result for small positive e is as follows. The non 
trivial zero of the beta function, G* ~ e/c, corresponds to an UV fixed point of 
the RG flow. The domain < G < G* is a weak coupling phase for gravity. At 
large distance G(p) — > and General Relativity is recovered. At small distances 
G(p) — > G* , and the fixed point governs the non-trivial short distance behavior of 
gravity at scales £ ~ Lpi anc k- The domain G > G* should correspond to a new 
strong coupling phase for gravity. Let me stress that this picture is valid for zero 
renormalized cosmological constant Ar. This should require the usual fine tuning of 
the bare cosmological constant in the regularized theory. 
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Fig. 19. The perturbative beta function for Einstein theory in D = 2 + e dimensions 

If this picture is valid for the four dimensional theory (e = 2), this would mean 
that gravity is "asymptotically safe" at short distance, thanks to the existence of a 
non-perturbative UV fixed point, which corresponds to a continuous phase transition 
between the weak coupling phase of gravity, with freely propagating gravitons, and a 
strong coupling phase, possibly with confined gravitons. 

However this interpretation of the perturbative calculation should be considered 
with some caution. The existence of a non-trivial UV fixed point which allows to 
construct a quantum field theory out of a perturbatively non-renormalizable theory 
is confirmed for the non-linear sigma models in 2 + e dimensions, while for gauge 
theories in 4 + e dimensions, a similar phenomenon has never been observed, and 
in numerical simulations the deconfinement transition (which separates the weak 
coupling deconfined phase from the strong coupling confined phase) is found to be 
first order in dimensions D > 5. The problem of the unboundness of the action, 
which do not raise inconsistencies in the perturbative renormalization calculations, 
still precludes a rigorous non-perturbative treatment of these questions. 

4-2. Simplicial 3D gravity 

One can use dynamical triangulations to discretize the functional integration over 
Euclidean 3-metrics (see sec. 2). Recently this approach has attracted much interest, 
in particular for numerical simulations. The basic idea is to start from some given 
initial 3d simplicial manifold, for instance with the topology of the 3-sphere S3, and 
to explore the space of combinatorially equivalent SM by performing the Alexander 
moves described in sec. 2, or rather the (1,4) and (2,3) moves (see Fig. 8), which are 
more efficient in numerical simulations. 

The possible discretized actions for the triangulation depend on the numbers of 
z-simplices Ni, i = 1,4, but using the geometric relation 4Ns — 2N 2 and the fact that 
the Euler characteristic for odd dimensional manifolds is zero 



x = N - iVi + N 2 - N 3 







(4.2.1) 
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one sees that there are only two quantities, for instance No (number of vertices) 
and Ns (number of tetraedra), which can vary independently. At variance with the 
two dimensional case, it is not possible to work in the canonical ensemble where the 
number if simplices is fixed. In fact, to go from a given small triangulation to another 
small one, one must pass in general through configurations with a large volume. 

Let me also mention that one can think of generalizations of matrix models to 
generate 3d triangulations. The simplest idea is to consider tensorial models (see 
Gross (1992), Ambj0rn et al. (1991)), where the field Tijk depends now of 3 indices, 
with a quartic interaction term (which is diagrammatically dual to a tetraedron) in 
the action 

^(T") — — Tijk Tkji g Tijk Tj n i T\ m k Ti mn . (4.2.2) 

However, the Feynman diagram for such models are in general not dual to 3d simplicial 
manifolds, but only to 3d pseudo-manifolds with non zero Euler characteristic. In 
addition, such models are not exactly solvable in the large N limit like the matrix 
models. 

In the DT approach the partition function is a sum over 3d triangulations T with 
a fixed topology, of the form 

Z(a,P) = J]e aJVo -' 9JV3 . (4.2.3) 

T 

This is the discretized form of the functional integral over 3d Euclidean metrics given 
by eq. (1.1.1-2). Indeed N3 corresponds to the volume of the manifold and No is 
linearly related to the total scalar curvature (an exact formula can be written using 
the deficit angle formula (2.4.3)). Large positive a favors positive total curvature, 
while large negative a favors negative total curvature. 

It is easy to show that < No < cst.N^. However, to show that this sum makes 
sense, the number of triangulations with a given topology and a fixed volume N3 must 
not grow faster than cst. N3 . In other words, the entropy of triangulations must be 
proportional to the volume. This ensures that the sum given by eq. (4.2.3) converges 
for f3 (the cosmological constant) large enough, and that in MC simulations one 
reaches an equilibrium configuration for simulation time t — > 00. 

Contrarily to the two dimensional case, there is no mathematical proof that such 
a bound holds! Some plausibility arguments have been given by Ambj0rn, Durhuus 
and Jonsson (Ambj0rn et al. (1991)). Numerical estimates (Agishtein and Migdal 
(1991b), Ambj0rn and Varsted (1991), Boulatov and Krzywicki (1991), Ambj0rn et 
al. (1992a)) for triangulations with the topology of the 3-sphere and up to N% ~ 10 4 
indicate that such a bound is valid and that Z(a,(3) defines a series (in e - * 9 ) with a 
finite radius of convergence. 

This result, which implies that the partition function Z(a,f3) exists for f3 > /3 c (a), 
and becomes singular at some critical value f3 c (a) of the bare cosmological constant, 
is a first non-trivial outcome of numerical simulations. As for the 2d case, one expects 
that at f3 c the average volume (N3) diverges and that the renormalized cosmological 
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constant vanishes. The next important issue is to understand this critical point, and 
whether it has some relationship with continuum 3d gravity. 

For this purpose, one must study how the observables of the model depend on a 
along the critical curve f3 = f3 c (a). The simplest observable is the total curvature 
(No). The numerical simulations show the existence of two different phases. For 
a > a C} where a c ~ 4.0, (No) is large and is found to be proportional to (N3) (cold 
phase). For a < a C} (No) is much smaller, and seems not to grow linearly with 
(N3), but rather as (No) oc (N%) s with 8 ~ .6 (hot phase). It is crucial to determine 
the nature of the transition between these two phases, in particular to see whether 
it is discontinuous (first order), or continuous, since in the latter case there should 
correspond a continuum field theory, with two relevant couplings (a and /3), in analogy 
with the results of the perturbative calculations (2 + e expansion) discussed above. 
Unfortunately, from the most recent MC simulations, an hysteresis cycle is observed 
around a c (see Fig. 20). This indicates that the transition is probably first order. 



< No/N 3 > 



.4 



hot phase 




cold phase 



3.0 3.5 4.0 4.5 5.0 OC 

Fig. 20. Hysteresis cycle for the curvature at the hot/cold phase transition for 3d DT 
(redrawn from Ambj0rn et al. (1992a), no error bars represented) 



Let me propose a tentative interpretation of these results. In the cold phase the 
curvature is large and positive. It is tempting, by analogy with the 2d case, to 
assume that in this case branched polymer configurations are dominant. Such BP 
configurations for 3d S^-like metric can be obtained very easily by starting from the 
smallest triangulation of S3 (the boundary of the 4d simplex), and by iterating the 
(1, 3) move, i.e. by expanding one tetraedron into three tetraedra, in randomly chosen 
tetraedra. We have no proof that such configurations maximize No , but it is easy to 
see that asymptotically No — §^3, which agrees with the numerical results for a > 4. 
(see fig. 20). 
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step 3 



Fig. 21. A singular 3-sphere 



In the hot phase the curvature is large and negative. We expect that in this case 
collapsed configurations such that a few links carry very large negative curvature 

1 12 

dominate. Let me give a simple example of such a configuration, such that No ~ iV 3 
(see fig. 21). Glue P tetraedra along a common edge to form a "pancake" (step 1). 
Then stack Q such pancakes to get a "ball" with an "equator" made of P links and 
with Q interior links between the "north" and "south" poles (step 2). Finally identify 
the north and south poles and the corresponding hemispheres to obtain a 3-sphere 
with Ns = PQ and No = P + Q (step 3). Such a configuration is in fact quite similar 
to the singular 2d triangulations that we discussed in subsec. 3.6. The two points 
with very large negative curvature are replaced here by two intertwinned loops. 

If this still speculative interpretation is correct, this would mean that the phase 
diagram for the 3d DT model is somehow similar to the phase diagram of 2d 
gravity in the strong coupling regime c ^> 1. For positive Newton's constant 
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the conformal instability would lead to branched polymer configurations with large 
positive curvature R oc a~ 2 (where a is some UV cut-off), while for negative Newton's 
constant collapsed configuration with large negative curvature localized along defects 
(points in 2d, loops in 3d) are dominant. It is not clear to me if there is space for a 
weak coupling phase where the curvature is small, and such that the Einstein theory 
is recovered in the classical limit. Of course we expect that numerical simulations 
with larger systems will lead to a better understanding of the hot and cold phases 
and of the nature of the transition. 

Finally let me mention that numerical simulations of Ising models on 3d dynamical 
triangulations have been performed recently (see Ambj0rn et al. (1992b)). The first 
results lead to qualitatively similar behavior when 3d gravity is coupled to that kind 
of matter fields. 

4-3. Simplicial 4D gravity 

The same techniques can be developed for Ad simplicial gravity. Ad simplicial 
manifolds are constructed by gluing 4-simplices along their 3d faces (tetraedra). 
Ergodic moves in the space of equivalent triangulations can be achieved by the 
three moves (1,5), (2,4) and (3,3). The condition that 4-simplices are glued along 
common 3-simplices implies 5N4 = 2N3. Using the fact that the boundary of 
the 4-ball formed by the 4-simplices which share a common vertex must be a 3- 
sphere (since the triangulation is a simplicial manifold), with Euler characteristic 
X = 0, one obtains, by summation over all the vertices, the additional constraint 
2N\ — 3N2 + AN3 — 5N4 = 0. These two constraints, together with the Euler 
characteristic identity \ = No — N\ + N2 — N3 + N4, imply that in a given topological 
class of triangulations only two quantities (for instance the number of vertices No and 
the volume N4) can vary independently. 

Numerical simulations have been performed by several groups for Ad DT with the 
topology of the 4-sphere (x = 2), with an action similar to that of eq. (4.2.3) 

S = -a No + /3iV 4 (4.3.1) 

for triangulations with up to N4 ~ 10 4 . The results are qualitatively similar to the 
3d case. I refer to Agishtein and Migdal (1992), and to Ambj0rn and Jurkiewicz 
(1992). One obtains a critical point at a finite value of /3, indicating that the entropy 
of Ad triangulations still grows linearly with the volume. As the curvature coupling 
a vary one goes from a hot phase, where (No) <C (^4), to a cold phase where (No) 
grows linearly with (N4). The simulations on large systems do not indicate hysteresis 
at the transition between the cold and hot phases. Therefore it is possible that the 
transition is continuous, at variance with the 3d case. One also observes that in 
the cold phase, for large a, (No) — .25 (N4), which is consistent with a branched 
polymer interpretation where No = N±/A. However, finite size effects are expected 
to larger in Ad (a lattice with 10 4 simplices has a diameter < 10 !) and it is clear 
that simulations on larger systems have to be performed, and will probably reveal 
unexpected phenomena. 
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4-4- 3D and 4D quantum Regge calculus 

I have no time to discuss here the results of numerical simulations using fixed 
triangulations with varying edges lengths. I refer to the recent reviews by Hamber 
(Hamber (1991a, b)). In these simulations also one observes a phase transition between 
a hot phase with large negative curvature and a cold phase with large positive 
curvature. The connection between these systems and the results of the DT is not 
yet very clear to me, and I think that the problems that I discussed for the two 
dimensional case still exist in higher dimensions. 



5. Non-perturbative problems in 2D quantum gravity 

In this last lecture I shall come back to two dimensional gravity and to matrix models, 
and discuss the issue of topology changes. This is interesting for two reasons. Firstly, 
topology changes are much simpler in two dimensions than in higher dimensions, 
since the topology is characterized by the genus, and explicit results can give hints 
about the higher dimensional models. Secondly, 2d gravity theories correspond to 
solutions of string theory in some simple backgrounds, and the topological expansion 
corresponds to the perturbative expansion around these backgrounds. Therefore one 
might hope to get a better understanding of string theories beyond perturbation 
theory through the matrix models. I shall mostly discuss here the one-matrix model, 
which allows to construct "pure gravity" (c = 0), and gravity coupled to non-unitary 
matter (described by (2, 5) minimal models in the classification of CFT). 

5.1. The continuum double scaling limit 

To fix notations, let me recall how the so-called double scaling limit (Brezin and 
Kazakov (1990), Douglas and Shenker (1990), Gross and Migdal (1990)) is obtained 
for the m = 2 critical point, which corresponds to pure gravity (see the lectures by 
E. Brezin in this volume and Gross et al. (1992) for details). The Hermitian matrix 
integral (3.3.1) 

d$ e- WTr (^)) (5.1.1) 
can be rewritten as an integral over the N eigenvalues (e.v.) 

r N 

Z = Jl[dX l A(X i ) 2 e- N E v (x i ) . A(Aj ) = JJ( Aj _ Aj ) . (5.1.2) 



1=1 Z<J 



In the planar limit N — > 00, this integral is dominated by the "master field" 
configuration corresponding to a continuous distribution of e.v's described by the 
local density po(X) } which extremizes the effective action (Brezin et al. (1978)) 

S eS [p] = [ dXp(X)V(X) - [ dX [dpp(X)p(p) In I A — // 1 (5.1.3) 
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subjected to the normalization constraints 



dAp(A) = 1 



p(\) > . 



(5.1.4) 



The first term in eq. (5.1.3) represents the external potential V, the second term 
corresponds to the repulsion between e.v's. The fluctuations of the e.v. density are 
subdominant (of order 1/iV 2 ), and the dominant contribution to the free energy is 
just 



F = ln(Z) ~ N 2 S eS [p 



(5.1.5) 



In the cases we shall be interested in, the potential V is such that the e.v's are localized 
along a single segment [a, b] on the real axis, and generically the e.v. density behaves 
at the end-points a and b as 



p(X) oc |A e - A] 1 / 2 



(5.1.6) 



(p is a deformation of the celebrated Wigner semi-circle distribution obtained for the 
quadratic potential V = A 2 /2). 

The m = 2 critical point is reached when the potential is such that the e.v. density 
becomes singular at one of the end-points, A e , and behaves as 



p(X) oc |A e - A| 3 / 2 . 



(5.1.7) 



For the cubic potential of eq. (3.3.2), V(A) = A 2 /2— gA 3 /3, this occurs at g c = 2-3 _3//4 . 
For g < g C} i.e. close to the critical point, p still has a square root singularity at the 
end-point \ e (g) < A c = A e (g c ), and the explicit solution shows that (A c — A e ) oc 

(g c -g) 1/2 . 




A, A, A, A, 

c c 

(a) g=g c (b) g<g c 

Fig. 22. The e.v. density at g c (a), and close to g c (b) 
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One obtains the double scaling limit by considering a small interval of width AX 
around the critical end-point A c , and by rescaling the parameters g ~ g c and N ~ oo so 
that there is only a finite average number of e.v's (n) ~ 7 -1 in this "critical interval". 
From eq. (5.1.7) this average number of e.v's in the critical interval [A c — AX, X c + AX] 
is easily estimated to be of order 

(n) ~ N dXp(X) ~ N (AX) 5 ' 2 . (5.1.8) 

J\ C -AX 

2 

Thus we require that N scales as AX ~ (iV 7 ) ' . The size of the critical region 
for the coupling constant g such that the end-point A e stays in the critical interval 

[A c - AX,X C + AX] is therefore \g - g c \ ~ Ag ~ (AX) 2 ~ (iV 7 )" 4/5 . 

Let us compare this with the continuum limit for the DT model (at fixed topology) 
discussed in subsec. (3.2). Since g is identified with the bare cosmological constant 
A of eq. (3.2.7), the short-distance cut-off a corresponds to the size of the critical 
region AX, and we can write the renormalized cosmological constant Ar, which will 
be denoted by to in the following, as 

t = Ar = a~ 2 (g c -g) , (5.1.9) 

while we rescale the e.v. position A as 

p = a' 1 (X - A c ) . (5.1.10) 

p corresponds to the "loop momentum", conjugate to the renormalized loop length 
£. To keep I/7 e.v's in the critical interval |p| < 1 amounts to rescale N (the total 
number of e.v's) as 



7 



x 5 / 2 N . (5.1.11) 



In the double scaling limit we let a — > while to, 7 (and p) are kept finite. From 
the KPZ scaling prediction (given by eq. (3.1.10) with 7^ = —1/2 for the c = pure 
gravity case), the term of order h of the topological expansion of the partition function 
F, which is given by the sum over connected closed surfaces with genus h, is expected 
to behave for g ~ g c as (t^ 2 j -2 )^ 1- ^ , and the topological expansion for F should 
become a series in the single scaling variable to 



00 



*X*o) = Y, J 2(h - 1} fh(to) ; fh(to) = (to) 5/2(1 h) h ■ (5.1.12) 



h=0 



Of course, in this simple case, the string coupling constant 7 can be reabsorbed in 
to- We keep it explicitly in order to simplify the discussion of the planar limit 7 — > 0. 
In this limit the number of e.v's in the critical interval becomes infinite, and their 
fluctuations can be neglected. 
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5.2. The Painleve I string equation 

It is known since some time that there is a deep connection between the statistics of 
e.v's of random matrices and some integrable systems (see Jimbo et al. (1980) and 
the book by Mehta (1991)). One way to understand this connection is to reformulate 
the integral over the e.v's as a problem of non-interacting fermions in one dimension. 
Introducing the orthonormal polynomials 7r n (A) for the measure dA e~ NV ( x \ 

deg(7r n ) = n , / dAe- JY ^ A ) 7 r m (A)7r n (A) = 6 mn , (5.2.1) 



the Vandermonde determinant is proportional to the Slater determinant for the first 
N polynomials 

A(Ai) oc det(vr J (A J )) , (5.2.2) 

and one can consider the "wave functions" 

tf„(A) = vr n (A) e-f v ^ = (A|n) (5.2.3) 

as one-particle states for a system of free fermions. The vacuum state |o) is obtained 
by filling the first N levels, since there are N e.v's. The correlations functions 
are encoded into the "position operator" Q } which acts by multiplication by A, 
7Tn(A) i— > A7r n (A), whose one-particle matrix elements are defined by 

Q nm = (n\Q\m) = f dAe- Jv ^ A )A 7 r n (A)7r m (A) . (5.2.4) 



Q is Hermitian and real and elementary properties of orthogonal polynomials show 
that only its first subdiagonals are non-zero, i.e. Q nm = if \n — m\ > 1. For instance 
the v.e.v. for the (Laplace transformed) one-loop operator W(\) = Tr((A — $) _1 ) is 
given by 

1 JV " 1 1 

(W(X)) = (o\j-q\o) = Y,( n \j-Q\ n ) ■ ( 5 - 2 - 5 ) 

The matrix elements of Q satisfy the so-called "recursion relations" (Itzykson and 
Zuber (1980), Chada et al. (1981)), which follows from the Heisenberg algebra 
satisfied by Q and the derivation operator P, n \— > -jp^ n . P is antisymmetric 
(P = P r ), and from eq. (5.2.3), P + fV(Q) is upper triangular. Therefore it can be 
written as 

P = ^(V'(Q) + -V'(Q).) (5.2.6) 

where ^4+(-) is the upper (lower) triangular part of the operator A ((m\A+\n) = 
(m\A\n) if m < n, otherwise). The commutation relation 



[P,Q] = i 



(5.2.7) 
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together with eq. (5.2.6), gives non-trivial equations for the matrix elements of 
Q } which allow to compute them recursively. This commutation relation can be 
considered as an isomonodromic deformation condition on the linear system obtained 
from eq. (5.2.6), which is of the form 

^ - = , (5.2.8) 

where \I/(A, n) is a 2-components field constructed out of the \I/ n 's (see Moore (1991)). 

In the double scaling limit, we may wonder which "states" \n) contribute to the e.v's 
density N p(X) = ir~ 1 lm( y (W(\ + «e))) for A in the critical interval | A — A c | < a, i.e. for 
|p| < 1. It appears that these are the states which are close to the "Fermi surface", and 
more precisely those states |m) with level number m such that N — m < [jV zlA) 1 . 
Therefore the number of such states, although small w.r.t. the total number N of 
states, is large w.r.t. the average number of fermions in the critical interval. In the 
double scaling limit the label m can be rescaled in a continuous parameter x 

-2 (-> m \ 



x = a-' (1--) . (5.2.9) 

For the c = critical point one can start from the cubic potential (deg(V) = 3). In 
the double scaling limit Q will become a Hermitian differential operator of degree 2 
(remember that only the first subdiagonals of Q are non-zero), and from eq. (5.2.6) 
P becomes an anti-Hermitian differential operator of at most degree 3. The explicit 
calculation shows that Q takes the form 

Q = 7 2 ^T " 2u(x,t ) (5.2.10) 
ox 

(Douglas (1990)) and the recursion relations for Q become the so-called "string 
equation" for it, which is the Painleve I equation 



7 2 d 2 u 
~6d~ 



x z 



+ u 1 = x + t . (5.2.11) 



The vacuum energy F is explicitly given by the integral (defined by a finite part 
prescription at oo) 



poo 

-2 



F(t ) = -7 _i / dxxu(x,t ) . (5.2.12) 
Jo 

Hence u depends only on x + to and — u is the second derivative of the partition 
function u(to) = —F"(to). The one-loop operator is (Banks et al. (1990)) 

(w(p)) = / dx(x| \x) . (5.2.13) 

Jo p-Q 
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In the "classical limit" 7 — s- 0, one should recover the results of the large N 
calculations. One has explicitly at leading order in 7 

u(x,t ) = \/x +t , ^0) = ~J~ 2 TP (to) 5/2 , 

15 

(Hp)) = -7" 1 I (tl /2 - p) [p + K 12 ) 1 ' 2 • (5-2.14) 



5.3. Non-perturbative properties of the string equation 

The Painleve I equation, that we rewrite by shifting x <— x + to as 

7 2 

— — u" -\- u = x , (5.3.1) 

is the simplest second-order differential equation whose solutions have the so-called 
Painleve property. Any solution of eq. (5.3.1) has an essential singularity at x = 00, 
and is a meromorphic function in the whole complex plane, with only double poles as 
"movable singularities" (i.e. singularities depending on the initial conditions). These 
double poles have the same residue, and in fact at any pole xo 

u(x) = ^— + O((x-x ) 2 ) . (5.3.2) 

(X — XQf 

The initial condition at +00, u(x) ~ y/x, and the Painleve equation are sufficient to 
specify uniquely all the terms of the asymptotic expansion of u at large x 



00 



u(x) = / y ttfc 7 x ~ 

k=0 

1 49 

iin = 1, ui = , uo = , ... (5.3.3) 

48' 4608' V ; 

However, this is not sufficient to fix the whole solution. Indeed, let us assume that 

there exists two solutions u and u with the same asymptotics at +00. Linearizing 
eq. (5.3.2) for the difference 8u = u — u we obtain at leading order 

- Su" + 2^Su ~ , (5.3.4) 
6 

whose solutions with growth slower than y/x are exponentially small 

Su(x) oc x" 1 / 8 e-^'t 2 ^ 6 ' 4 . (5.3.5) 
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Fig. 23. Schematic view of a real solution of the Pj string equation 



In fact a closer analysis of the Painleve I equation shows that it has an infinite 
family of solutions with the same asymptotic expansion at +00 (see Boutroux (1913) 
and the book by Hille (1976)). These solutions differ, for instance, by the position 
of the largest double pole on the real axis. The difference between two solutions 
is exponentially small, from eq. (5.3.5), and is invisible in perturbation theory. This 
kind of phenomenon is typical of non-Borel summable series (see below). The analysis 
of the asymptotic behavior of u for negative x shows that the real solutions must 
necessarily have an infinite series of double poles along the negative real axis. Indeed, 
it is natural to rescale 

X = -x 5/4 , U(X) = x~ 1/2 u(x) , (5.3.6) 
5 

so that the Painleve I equation becomes an equation for U(X) } which is asymptotic 
for \X\ — > 00 to an elliptic equation (I have set 7 = 1) 

U" ~ QU 2 + 6 = - V - + = OiX- 1 ) . (5.3.7) 

The generic solution when the r.h.s vanishes is a Weierstrass function p{X) with a 
doubly periodic lattice of double poles in the complex plane. A careful analysis shows 
that if u is asymptotic to iy/x as x — > +00, generically this analytic asymptotics 
is valid in two-fifth of the complex plane, namely for —2tt/5 < Arg(x) < 2tt/5. In 
the remaining three sectors, which include the negative real axis, u has necessarily 
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an infinite number of poles. There is a unique solution, called the "triply truncated 
solution", which is asymptotic to iy/x in the larger "sector of analyticity" —4tt/5 < 
Arg(x) < 4tt/5. This solution is real, but still has poles on the negative real axis, and 
has the wrong sign asymptotics u ~ — \fx for our problem. 

Going back to the "physical solutions" of the string equation, which must satisfy 
u(x) ~ y/x, we would like to understand the origin of these non-perturbative terms, 
and to know which physical constraint, if any, allows to choose between the different 
real solutions of eq. (5.3.1). The unavoidable existence of an infinite number of double 
poles of u for negative x is a serious problem for defining the resolvent (x\(p — |y), 
and the continuum limit for the loop operators (David (1990)). 

5.4- Divergent series and Borel summability 

Asymptotic but divergent series expansions are very common in quantum mechanics 
and in quantum field theories. A mathematical recipe to make sense out of these series 
is the Borel summation method. I refer to Zinn- Justin (1983) and to Le Guillou and 
Zinn- Justin (1990) for details. It relies on the following Theorem: 

If the function f{x) is analytic in the sector S = {x : |Arg(x)| < a7r/2, \x\ < est.}, 
and has a Taylor expansion f(x) = ^^_ fk% k , such that the following bound on 
the rest of the expansion holds uniformly for any x in S and at any order N , for 
some f3 < a 



N 



/(*) " Y,f kxk ^ BA- N (/3N) 



\x N+1 



(5.4.1) 



A;=0 



then f is unique. It can be defined from its Borel transform (of order f3) 



OO j, 

which is analytic in the domain {b : |6| < A or |Arg(6)| < a/2} by the Borel 
summation formula 



[•OO 

f{x) = dt e"* B{xtP) . (5.4.3) 

Jo 

In fact, at least for f3 = I, it is sufficient to have a sector with a = f3 (Nevanlinna- 
Sokal theorem). 

In general, in quantum mechanics and in quantum field theories, one has f3 = 1 
(but f3 may be smaller for theories with fermions, such as QED). The singularity 
of the Borel transform closest to the origin, bo, is usually real and negative, and 
from eq. (5.4.3) the terms of the series fk grow factorially, but with alternate signs, 
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as fk oc kl( — l) k \bo\~ k . These singularities are in fact related to the existence 
of possibly complex Euclidean solutions of the equation of motion with finite 
action, called instantons. Additional singularities of the Borel transform, usually 
called renormalons, arise from short distance divergences (for strictly renormalizable 
theories), and from large distance divergences (for perturbatively massless theories). 
I shall not discuss these renormalon singularities here (see Le Guillou and Zinn- Justin 
(1990)). As a crude argument, let us consider the Euclidean functional integral over 
field configurations <f> for a theory with action S(<f>), and coupling constant g (which 
plays the role of Ti). The partition function can be rewritten as 

Z(g) = / d[<f>] e~ s ^/3 = / e -b/g B ^ ^ ( 5A4 ) 



■ioo 2i7r 



Wl 



ith 



e m = J d ^sw^ ■ (5A5) 

The function B is expected be analytic for negative 6, but as Re(6) increases, it will 
become singular for b = bo = S(<f>o), where cf>o is the extremal configuration with lowest 
action. By appropriate contour deformations, one can analytically extend B into a 
function of b with cuts starting at some h{ = S(<f>i) corresponding to saddle points 
configurations (pi of S, with finite action. These classical solutions are the instantons. 
In the semi-classical limit g — > 0, the perturbative expansion for Z(g) is obtained 
by moving to the right the contour over b in eq. (5.4.5). The dominant contribution 
is given by the first cut, associated to the lowest action configuration — the classical 
vacuum — with action So • The Borel transform of the perturbative expansion around 
the vacuum is simply the discontinuity of B around the first cut 

1 f°° 
Bo(b) = - lm(B(b)) , Z(g) = / db e~ b ^ B (b) . (5.4.6) 

However, if the <f> integration for B "catches" other instantons, the contributions of the 
other cuts associated to these instantons have to be added in eq. (5.4.6). Through the 
Stokes phenomenon, they will contribute to cuts in the b plane for the perturbative 
Borel transform Bo(b). 

For simple theories, such as the massive scalar \<f> 4 theory (with A > and in 
dimensions D < 4), there are no real instantons with positive action, the perturbative 
Borel transform Bo is well-defined for b > 0, and the theory is Borel summable. The 
closest singularity of Bo(b) is negative. The corresponding instanton is an imaginary 
solution of the equations of motion, or equivalently a real solution for the upside-down 
potential (A < 0). However in theories with inequivalent topological sectors, such as 
non-abelian gauge theories, real instantons which describe tunneling between these 
vacua do contribute to the functional integral, and have to be taken into account. 
They will give non-perturbative contributions of order e~ Sinst / g . 
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Another situation where instantons are crucial is of course the case of theories with 
an action unbounded from below. Instantons describing the decay (through quantum 
tunneling) of the unstable perturbative vacuum give singularities of Bo(b) on the 
positive real axis, but in this case there is in general no way to sum consistently 
perturbation theory. A simple O-dimensional example is given by the integral 

Z(g) = I dA *-*lWl> = f" gt W (5A7) 

It describes a single particle in a cubic potential. The perturbative expansion is 
obtained by expanding around the metastable vacuum Ao = 0, but there is a second 
saddle point at Ai = 1/ ^fg corresponding to the particle at the top of the wall. This 
O-dimensional instanton has an action S\ = l/(6g), it gives a singularity of the Borel 
transform at b = S\, and the coefficient of order g k behaves as k\6 k (but with no 
alternate signs), so that the series of eq. (5.4.7) is not Borel summable. This is of 
course a consequence of the divergence of the integral as A — > +oo. 



s 




Fig. 24. A zero-dimensional instanton 

The above integral defines nothing but the number of Feynman diagrams for a 
bosonic theory with cubic interaction. The fact that the total number of Feynman 
diagram with k internal loops grows like k\ (constant)* is a general feature of local 
quantum field theories. Moreover (see Rivasseau (1991b)), it is possible to show 
that the amplitude of a "typical" diagram at order k grows also as (constant^ (not 
taking into account the UV and IR singularities that occur in some specific subclasses 
of diagrams, and which give rise to the renormalon singularities mentioned above). 
Thus it is the instanton with smallest action which governs the large order behavior 
of each term of the perturbative expansion. 

5.5. Non-perturbative effects in 2D gravity and string theories 

For the c = pure gravity theory, the form of the non-perturbative terms, given 
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by eq. (5.3.5), indicates that the term of order k of the topological expansion will 
not grow like k\, but like (2k)\ Indeed, the expansion parameter is 7 2 , but the non- 
perturbative terms behave as exp(— 7 -1 est.), not as exp(— 7~ 2 est.). This suggests to 
use a f3 = 2 Borel transform w.r.t. the string coupling constant 7 2 , that is an ordinary 
Borel transform (/3 = 1) w.r.t. 7. The free energy F(t) is for instance 



00 

1 .5/4 



F(t) ~ / dbe~^~ f B(b) . (5.5.1) 



The existence of a one-parameter family of solutions of the string equation, with the 
same perturbative expansion, but which differ (from eq. (5.3.5)) by terms of order 
8F oc t -5//8 exp ( — 7 _1 6it 5//4 ), with b\ = 8y/3/5, implies that B(b) has a singularity 
at 61, and behaves as (b\ — 6) _1//2 , so that the perturbative expansion (5.1.12) for F 
is not Borel summable, and that each term fk has the same sign and grows like 

/* oc T(2k - |) b~ 2k . (5.5.2) 

It is possible to understand this (2k)\ large order behavior, and the corresponding 
singularity of the Borel transform at 61, in terms of an instanton effect, but in this 
case only one among the N eigenvalues is tunneling (David (1991)). Indeed, in order 
to reach the m = 2 critical point, one must start from a potential V unbounded 
from below in eq. (5.1.2), and for N large but finite the integral over the N real 
eigenvalues is divergent (but it can be used as a formal generating functional of 
the topological expansion). Looking at the variation of the effective action given 
by eq. (5.1.3) when one eigenvalue is moved from its equilibrium position to some 
A outside the support of e.v's, (this amounts to change the distribution density by 
Sp(fi) = iV -1 (<*)(/i — A) — 8(p — Ao) with Ao G Supp(/?)), one obtains the effective 
potential T(A) for one eigenvalue in the background created by the N — 1 remaining 
e.v's, which are supposed to stay at equilibrium. Y is given by 



r(A) = n 



V(X) - 2 / d/i / o(/i) In I A - /i| 



(5.5.3) 



The saddle point equation for p is equivalent to the equilibrium condition that the 
force 



/(A) = -T'(A) = N 



V'(X) + 2 J- d/i 



Pit*) 
X — /j, 



(5.5.4) 



must vanish wherever p(X) ^ 0. This implies that T = constant on the support 
of p. This condition, together with the normalization and positivity constraints of 
eq. (5.1.4), fixes the eigenvalue distribution p. 

In the double scaling limit, only the critical neighborhood of the e.v. distribution 
end-point is important. If we perform the rescalings given by eqs. (5.1.9-11) and 
take 7 — )• 0, i.e. the number of e.v's in the critical region to be large, we obtain the 
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force, which is given by f(p) = 2He(w(p)) , with the loop operator (w(p)) given by 
eq. (2.2.13). Integrating w.r.t. p we obtain the expression for the effective potential 

T(p) = Re (V 1 ± (3tJ /2 - p) (p + 2tJ /2 ) 3/2 ) . (5.5.5) 

The shape of T, together with that of the density pip) = Im(tz;(p — ie))/27r, is depicted 
for t > on Fig. 25. 




Fig. 25. One e.v. effective potential T (black) and the density p (dashed). 

Perturbation theory is constructed by taking into account the fluctuations of the 
e.v's around the "classical vacuum" where all the e.v's are on the constant potential 

1/2 

half-line p < po = —2t One may consider the instanton configuration where one 

1/2. 

eigenvalue sits at the top of the potential well p\ = t , while the others stay in the 
classical vacuum. At leading order in 7 one may neglect the backreaction of this e.v. 
on the background configuration, and the action for this configuration is 

r inst = r( Pl ) = t" 1 ^/ 4 , (5-5.6) 

which corresponds exactly to the strength of the non-perturbative effects given by 
eq. (5.3.5) in the string equation. 

In order to estimate the corrections to this behavior, for instance the t 5//8 power 
correction, one needs to take into account the backreaction of the e.v. on the bulk. It is 
also expected that higher orders non-perturbative terms have a similar interpretation 
in terms of "multi-instantons" , i.e. of several e.v's sitting on the top of the effective 
potential. 
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This (2k)\ large order behavior is expected to be a very general feature of string 
theories (Shenker (1991)). Indeed, in the same way as the behavior of the single 
integrals corresponding to O-dimensional field theories (such as that of eq. (5.4.7)) 
allow to recover the k\ growth of the number of Feynman diagrams with k loops 
in quantum field theories, simple matrix models allow to estimate the volume of the 
moduli space Ai g of closed Riemann surfaces with genus g. For instance, pure gravity 
is equivalent to the Liouville+ghosts system and the term of order k of the topological 
expansion of F can be written as 

F k = / d/i[m] ^LiouvilleM ZghostsM , (5.5.7) 

JM k 

where djj\m\ is the natural measure (the Weil-Peterson measure) over the moduli 
m, and Z...[m] are the functional determinants of the fields living on the surface. 
For the c < 1 theories this integral is expected not to diverge on the boundary of 
Mk (but tachyonic divergences will occur for c > 1), and to be of the order of 
Vol(Aik) = Xmj. dfi[m]. This volume can be estimated directly by some random 
matrix models that produce triangulations of the moduli spaces (Harer and Zagier 
(1886), Penner (1987)), and by topological field theory (Witten (1990)). 

It is expected that for the other c < 1 non-critical string models, the non- 
perturbative effects can also be understood in terms of one e.v. instantons. This has 
been checked explicitly for the c = 1 model. This means that the field configurations 
that govern the non-perturbative effects in string theory cannot be understood in 
terms of the low energy field-like modes of the string (tachyon, graviton, dilaton, . . .), 
which would give a k\ behavior, but are really "stringy" collective excitations. 

5.6. Various stabilization proposals 

We have seen on the example of the m = 2 critical point of the one-matrix model 
that instantons correspond to the tunneling of e.v's into the infinitely deep potential 
well. This means that the "classical vacua" is metastable, and that the theory has 
no real stable ground state. This observation can be extended to the matrix models 
which describe 2d gravity coupled to c < 1 unitary matter (Ginsparg and Zinn- Justin 
(1991). Some non-unitary models appear to be stable, and the topological expansion 
still behaves as (2k)\ } but is Borel summable, because of sign alternations of the 
matter partition function Z matter with the genus k. Thus, for the physically relevant 
models of non-critical strings, non-perturbative effects point out a serious problem 
of "unboundness of the action", not a simple tunneling between inequivalent — but 
stable — vacua as for gauge theories. 

Several attempts to modify the formulation of these theories, without changing 
their perturbative expansion, have been made, in order to "stabilize" the functional 
integral. I shall discuss briefly some of them. 

Stochastic quantization 

Marinari and Parisi (1990) suggested to treat the one-matrix model as a one 
dimensional supersymmetric matrix model, through stochastic quantization. It is 
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known that in general, the average of any quantity Q($) w.r.t. a probability density 
of the form d$ exp( — V($)) can be rewritten as the ground state expectation value 



of the observable Q($) for the quantum mechanical system with Hamiltonian 



with P = id/d<&. The time variable for this quantum mechanical system is the 
Langevin time when one uses stochastic quantization to compute (Q) } and Hb is the 
associated Fokker-Planck Hamiltonian. It is the restriction to the bosonic sector of 
a supersymmetric Hamiltonian H = Q 2 , where Q is the SUSY generator. When 
the potential is bounded from below, so that the probability density in eq. (5.6.1) 
exists, SUSY is unbroken, the ground state wave function is \I/($) = exp( — V($)/2) 
and has zero energy. If V in unbounded from below, the Hamiltonian Hb still makes 
sense, and it was proposed by Greensite and Halpern (1984) that it allows to define 
a quantum theory for systems with classical action unbounded from below. When 
applied to the one-matrix model, one obtains indeed a stable one-dimensional matrix 
model with variable $(t), where t is the Langevin time. It has the same topological 
expansion than the initial unstable zero-dimensional matrix model, but differs by 
non-perturbative contributions. These corrections come from the fact that, when the 
potential is unbounded from below, SUSY is spontaneously broken, and the vacuum 
energy is positive, but exponentially small. 

The properties of this model have been studied by several groups. I refer to Karliner 
and Migdal (1990), Ambj0rn et al. (1990), Miramontes et al. (1991) and Dabholkar 
(1990) for details. The main question raised by this proposal is that the KdV flow 
structure, the loop equations and the Virasoro constraints of the original model have 
disappeared. It is not clear what replaces them, and what is the interpretation of this 
new theory in terms of strings. 

Complex matrix models. 

In Dalley et al. (1992) it is shown that by using complex (non-Hermitian) matrices 
M and an action of the form Tr (F(MM')) one can construct well-defined models 
which have, in the double scaling limit, the same perturbative expansion than that of 
the ordinary model. The basic reason for this phenomenon is that, when expressed in 
terms of an integral over the eigenvalues of the matrix MM^ , the partition function 
can be written exactly as that of eq. (5.1.2), with the restriction that the e.v's must 
be positive (see Morris (1991)). This allows, when adjusting the potential V, to 
put a "wall" at a finite distance from the critical end-point of the e.v. distribution, 
which forbids the e.v's to fall, but which is invisible in perturbation theory. The 
string susceptibility u(x) is now found to obey a string equation which is related 
to the KdV hierarchy, but which is now 4 th order instead of 2 nd order. A problem 



(Q) 




(5.6.1) 





(5.6.2) 
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with this approach is that an additional non-perturbative parameter appears, which 
corresponds to the distance between the wall and the end-point. The loop equations 
are also broken by a term related to the disappearance of one of the Virasoro 
constraints. 

Complex solutions 

Finally, let me mention that one can stay with the one-matrix model and define the 
e.v. integral given by eq. (5.1.2) by choosing as integration contour not the real 
axis but a complex contour such that as |A| — > oo the potential V(A) is real and 
positive. For the cubic potential this amounts to take, for instance, the contour C+ 
going from — oo to e 17r / 3 oo. See David (1990), Sylvestrov and Yelkovsky (1990), David 
(1991). It is possible to show that with this choice, in the double scaling limit, one 
recovers a triply truncated complex solution of the Painleve I string equation, which 
is asymptotics to \fx for —6tt/5 < Arg(x) < 27r/5, and which has an infinite network 
of double poles in the remaining sector 2tt/5 < Arg(x) < 4tt/5. This solution, and 
the solution associated to the complex conjugate contour C_, can be obtained from 
the perturbative series by a well-defined Borel summation prescription (eq. (5.5.1)), 
where the contour of integration over the Borel variable b is taken above (or under) the 
instanton singularities located on the real axis. Of course this summation prescription 
violates the physical prescription of reality of the amplitudes. 

6. Conclusion 

I have not discussed many important issues in these lectures. The c = 1 matrix 
models solutions, and their relationship with 2d string theories and other conformal 
field theories, for instance those describing 2d black holes, have raised much interest 
during the last two years: see for instance Gross et al. (1992) and Polyakov (1992). 
Another important problem is to extend what has been done for 2d gravity to 2d 
supergravity theories, with the aim of obtaining exact solutions and non-perturbative 
information for superstring theories. Some progresses have been made recently (see 
the lectures by L. Alvarez-Gaume), but one does not have yet a simple geometrical 
discretization procedure of superspace similar to random lattices for ordinary space. 
Another interesting problem would be to construct discrete models and/or exact 
solutions corresponding to space-time with Lorentzian signature ( — |-)- 

For higher dimensions, in my opinion much remains to be done. I hope that in the 
next years large scale numerical simulations and new theoretical developments will 
help us to understand the nature of the different phases of these higher dimensional 
models, the importance of conformal modes, the connection between the different 
lattice formulations, and with the topological theories in three and higher dimensions. 
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